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SUMMARY 

This  report  contains  an  analysis  of  simultaneous  ground-based  measurements  of  infrared  (3-5/uw  and  8-12 nm)  and 
visible  (photopic)  atmospheric  extinction  The  measurements  were  taken  at  hourly  intervals  over  three  3-month  periods  in 
the  Netherlands,  as  part  of  the  NATO  OPAQUE  program.  The  measurements  were  made  by  the  Physics  Group  of  the 
Physics  Labora'ory  of  the  National  Defense  Research  Organization-TNO.  at  Ypenberg  Air  Base,  Netherlands.  The  analysis 
is  based  on  preliminary  data. 

A  scheme  for  sorting  the  infrared  aerosol  particle  extinctions  was  devised,  taking  into  consideration  the  form  of  the 
observed  parameter  relationships,  the  expected  measurement  uncertainty,  and  the  theoretically  predicted  conditions  associ¬ 
ated  with  high  infrared  aerosol  extinction.  The  data  were  sorted  into  an  upper  bin,  consisting  primarily  of  mist  to  fog 
cases,  and  a  lower  bin,  consisting  primarily  of  clear  to  haze  cases.  Here  "mist"  is  defined  as  a  suspension  of  water  droplets 
or  wet  hygroscopic  particles  (not  a  drizzle)  causing  reduced  visibility  greater  than  1  km.  When  the  visibility  is  reduced  in 
this  manner  to  I  km  or  less,  the  condition  is  defined  to  be  fog. 

The  upper  (mist)  bin,  defined  as  those  points  with  visual  extinction  coefficient  greater  than  1  km  1  (visibility  less 
than  3  km)  and  relative  humidity  of  94%  or  greater,  included  essentially  all  of  the  high  infrared  aerosol  extinction  values. 
Using  both  visual  extinction  and  relative  humidity  thresholds,  rather  than  either  threshold  alone,  yielded  much  more 
effective  sorting  of  the  infrared  extinctions.  These  thresholds  are  based  on  the  preliminary  data  release.  Re'ative  humidi¬ 
ties  near  94%  have  been  corrected  to  near  99%  in  edited  versions  of  the  data  base,  released  after  this  report  was  written. 

In  the  lower  bin  (Iwe)  category,  the  infrared  aerosol  extinction  (due  to  aerosol  particle  scattering  and  absorption) 
was  consistently  low  relative  to  the  extinction  caused  by  air  molecule  scattering  and  absorption  and  water  vapor  absorption. 
The  aerosol  extinction  in  the  3-5 fj.m  and  8-12^/«  regions  averaged  near  .02  km'1  in  the  summer  and  .01  km~!  in  the 
winter  These  values  compare  reasonabiy  with  LQWTKAN  model  estimates.  Thus  for  these  conditions,  i.e.  visual  extinc¬ 
tion  of  1  km  1  or  less  or  relative  humidity  less  than  94%.  the  aerosol  extinction  was  consistently  low. 

In  the  upper  bin  (mist)  category,  the  infrared  aerosol  particle  extinction  was  found  to  vary  over  a  wide  range,  with  a 
median  extinction  of  about  .1  km~!  and  with  a  significant  number  of  extinctions  of  1  to  10  km'1  or  higher.  The  aerosol 
extinctions  exceeded  the  air  molecule  and  water  vapor  extinction  only  about  20%  of  the  time  and  exceeded  1  km'1  only 
about  I0n/  of  the  nme  These  high  extinctions,  occurring  10-20%  of  the  lime,  compare  well  with  the  LOWTRAN  fog 
models  The  incidence  of  the  high  extinctions  wax  distributed  over  the  whole  range  of  visual  extinction  in  the  mist  bin. 
and  persistence  calculations  indicate  that  the  infrared  aerosol  extinctions  have  low  persistence  relative  to  the  visual  extinc¬ 
tion  and  relative  to  the  duration  of  the  mist  episodes. 

Thus  for  this  upper  bin  (mist)  category,  the  infrared  aer.iso*  extinction  reaches  significantly  high  magnitudes  less 
often  than  expected,  and  persis's  only  over  short  time  intervals.  Measurements  in  addition  to  visible  extinction  coefficient 
and  relative  humidity  would  be  required  to  predict  the  incidence  of  high  infrared  extinction  within  mist  episodes.  !n  the 
absence  o*'  this  information  one  can  still  make  statistical  estimates  of  occurrence  probabilities  using  cumulative  frequency 
distributions  of  the  infrared  to  visible  aerosol  extinction  ratio 
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An  Analysis  of  Infrared  and  Visible 
Atmospheric  Extinction  Measurements  in  Europe 


Janet  E.  Shields 


1.  INTRODUCTION 

As  electro-optical  systems  are  used  for  a  growing 
variety  of  applications,  there  is  an  increasing  requirement 
for  understanding  and  predicting  the  atmospheric 
influences  on  system  performance  In  support  of  ms 
requirement,  the  Visibility  Laboratory  has  maintained  up. 
extensive  program  of  visible  optical  and  meteorolog.-ai 
measurements  and  has  been  involved  in  the  analysis  of 
data  taken  under  the  NATO  p.ugram  OPAQUE  (Optical 
Atmospheric  Quantities  in  Europe!,  documented  in  Fenn 
(1978)  The  OPAQUE  program  consisted  primarily  of  a 
series  of  ground-based  optical  measurements  taken  in  the 
visible  and  infrared  portions  of  the  spectrum,  along  with 
approximately  s.mul'aneous  meteorological  measurements. 
These  measurements  were  recoiJed  ut  several  stations 
throughout  Europe  at  hourly  intervals  over  a  two  year 
period.  This  report  deals  with  an  analysis  of  some  of  the 
infrared  data  taken  during  Project  OPAQUE. 

There  are  a  variety  of  optical  infrared  sensors  in  use 
today,,  whose  performance  is  partly  dependent  on  the 
atmospheric  infrared  extinction  Mie  calculations  based  on 
model  particle  sice  distributions  result  in  the  following 
trends  [Shettle  and  Fenn  (1979)]  When  ihe  atmospheric 
conditions  range  from  clear  to  hazy,  aerosol  extinctions 
are  generally  smaller  at  infrared  wavelengths  than  at  visi¬ 
ble  wavelengths  Except  in  very  clear  air,  this  will  result 
in  better  optical  propagation  within  the  infrared  "window" 
regions  (about  3-5/xw  and  8-13/zm)  than  at  visible 
wavelengths.  However  in  the  presence  of  large  droplets  of 
size  about  equal  to  the  wavelength,  such  as  may  occur  in 
mist  and  fog,  the  infrared  aerosol  extinction  can  become 
greater  than  the  visible  aerosol  extinction,  resulting  in 
poorer  optical  propagation  of  infrared  wavelengths 

Unlike  the  visible  extinction,  the  infrared  extinction 
is  dnTicuii  and  expensive  to  measure.  Additionally,  the 
infrared  extinction  is  difficult  io  estimate  from  visible 
measurements,  since  the  infrared  to  visible  extinction  rela¬ 
tionship  is  influenced  by  the  aerosol  particle  size  distribu¬ 
tion  and  index  of  refraction,  which  are  also  difficult  to 
measure  As  a  result,  there  is  interest  in  evaluating  how 
well  present  models  estimate  the  infrared  extinction  and  in 
developing  improved  empirical  means  of  estimating  the 
infrared  extinction 


1.1  Goals  of  the  Analysis 

This  analysis  of  the  OPAQUE  data  was  directed 
toward  two  goals.  The  first  goal  was  to  evaluate  the  varia¬ 
tions  in  the  measured  data  relative  to  model  predictions. 
In  general,  haziness  and  the  associated  optical  changes, 
such  as  a  decrease  in  the  radiance  transmittance,  are 
caused  by  increases  in  the  aerosol  particle  content  of  the 
atmosphere.  This  aerosol  includes  droplets  which  may  be 
composed  of  varying  fractions  of  water,  sea  salt,  soot,  and 
so  on,  as  well  as  solid  particles  such  as  dust.  Measure¬ 
ments  have  been  made  of  the  size  distribution  of  these 
particles,  their  chemical  content,  and  their  indices  of 
refraction.  Most  models  are  based  on  Mie  calculations 
whose  inputs  are  the  most  likely  aerosol  particle  size  dis¬ 
tributions  and  the  most  likely  refractive  index  distributions 
as  a  function  of  wavelength.  Since  measured  aerosol 
characteristics  in  the  atmosphere  vary,  it  is  to  be  expected 
that  the  transmittance  also  should  vary.  The  extensive 
OPAQUE  data  set,  which  clearly  illustrates  these  typical 
radiance  transmittance  variations,  should  be  valuable  for 
evaluating  how  well  the  models  predict  them. 

The  second  goal  of  the  analysis  was  to  develop 
improved  techniques  for  estimating  the  infrared  transmit¬ 
tance  based  on  more  conventional  measurements  and 
observations.  Like  the  infrared  transmittance,  the  visible 
transmittance  is  also  influenced  by  the  aerosols.  The 
meteorological  quantities  are  also  related  to  the  aerosols. 
For  example,  aerosol  particle  size  is  influenced  by  relative 
humidity,,  and  aerosol  content  is  influenced  by  airmass  his¬ 
tory  Since  the  visible,  meteorological,  and  infrared 
regimes  affect  or  are  affected  by  the  aerosols,  it  might  be 
possible  to  develop  an  adequate  empirical  model  to  relate 
the  infrared  radiance  transmittance  directly  to  the  visible 
and  meteorological  quantities. 

2.  THEORY  AND  BACKGROUND 

2.1  Definition  of  Terms 

In  this  analysis,  the  parameter  of  interest  is  the  radi¬ 
ance  transmittance  in  the  infrared.  The  trarsmittance, 
which  is  a  function  of  path  length,  is  defined  at 
wavelength  a  as  the  proportion  cf  the  incident  beam 
which  remains  after  traveling  over  a  given  path. 


1 


rK(r)  -  /V,(r)//V.( 0) 


(2-1) 


where  Vk(r)  and  A\(0)  are  the  radiances  at  wavelength  A 
and  path  lengths  r  and  0.  Equation  (2-1)  is  correct  for 
monochromatic  radiation.  Broadband  sensors,  such  as 
those  utilized  during  Project  OPAQUE,  measure  a 
weighted  transmittance, 
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where  Nk  is  the  radiance  of  the  source  (in  this  case  a  650° 
blackbody)  and  Rk  is  the  spectral  response  of  the  sensor. 

The  monochromatic  transmittance  at  wavelength  A 
may  also  be  defined  in  terms  of  the  monochromatic 
extinction  coefficient  a,, 


vapor.  TMoi  is  the  transmittance  due  to  the  air  molecule 
components,  such  as  carbon  dioxide.  It  is  a  function  of 
temperature.  (For  convenience,  these  transmittances  and 
the  associated  extinction  coefficients  will  be  designated  by 
the  adjectives  "aerosol",  "water  vapor",  and  "molecular".) 
All  of  these  components  include  both  scattering  and 
absorption.  Although  Eq.  (2-6)  is  strictly  correct  only  for 
monochromatic  radiance,  treating  each  of  the  components 
as  an  integral  analogous  to  Eq.  (2-2)  results  in  errors  to 
the  total  transmittance  which  arc  usually  less  than  1%  in 
the  window  regions.  [Shettle  (1978a)]. 

Most  of  the  uncertainty  in  the  estimation  of  infrared 
transmittance  from  model  inputs  is  due  to  the  uncertainty 
in  estimating  the  aerosol  attenuation.  The  general  con¬ 
sensus  [Nilsson  (1979)  anti  Roberts  (1976)]  seems  to  be 
that  given  measurements  of  temperature  and  humidity, 
the  water  vapor  and  molecular  extinctions  are  well  handled 
by  computer  programs,  such  as  LOWTRAN  5  [Kncizys 
dal.  (1980)].  This  analysis  deals  with  the  aerosol 
transmittanre  and  the  associated  aerosol  extinction. 
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2.2  Related  Theoretical  Studies 


The  effective  broadband  extinction  coefficient  is 
defined  by  either  of  two  non-equivalent  definitions 
depending  on  the  application.  For  direct  broadband  meas¬ 
urements  of  absorption  and  scattering,  the  instrument 
response  depends  on  an  effective  extinction  coefficient 
defined  by 
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However  for  direct  broadband  measurements  of  transmit¬ 
tance  or  transmitted  light,  the  instrument  response 
depends  on  an  effective  extinction  coefficient  defined  by 
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(with  T  defined  in  Eq  2-2),  which  may  be  range  depen¬ 
dent.  This  range  dependence  occurs  when  the  attenuation 
is  wavelength  selective  within  the  band,  because  the  spec¬ 
tral  content  of  the  light  is  in  that  case  range  dependent. 
In  the  case  of  aerosol  extinction,  the  range  dependence  is 
very  slight,  and  Eqs.  (2-4)  and  (2-5)  are  nearly  equivalent. 

The  transmittance  may  be  represented  as  the  pro¬ 
duct  of  three  components, 
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Ta„  is  the  transmittance  due  to  the  aerosol  particles,  i.c. 
the  wet  and  dry  particles  in  the  atmosphere.  is  the 

transmittance  due  to  molecular  water  vapor.  It  is  a  func¬ 
tion  of  temperature  and  dewpoint  temperature  or  water 


There  has  been  a  great  deal  of  work  on  the  various 
theoretical  aspects  of  predicting  aerosol  particle  extinction. 
Two  rather  well  thought  out  works  are  those  of  Shettle 
and  Fenn  (1979)  and  Nilsson  (1979),  These  papers  dis¬ 
cuss  the  most  probable  aerosol  characteristics  such  as  par¬ 
ticle  size,  and  show  the  resulting  computed  extinction 
characteristics.  Pinnick  d  at.  (1979)  report  similar  calcu¬ 
lations  based  on  the  measured  drop  size  properties  of  fogs. 

In  general,  the  results  of  these  calculations  may  be 
summarized  as  follows.  If  the  relative  particle  size  distri¬ 
bution  and  refractive  index  distribution  ot  the  aerosol 
remain  constant,  and  only  the  number  density  changes, 
then  the  magnitude  of  the  aerosol  extinction  coefficient 
wiil  vary,  but  the  spectral  relationships  will  be  fixed.  That 
is,  the  ratio  of  extinction  coefficient  at  two  given 
wavelengths  will  be  a  constant.  If  relative  particle  size  dis¬ 
tribution  or  refractive  index  distribution  are  allowed  to 
change,  then  the  spectral  relationships  will  change. 

The  visible  extinction  should  be  more  strongly 
affected  by  the  sub-micron  region  of  the  particle  size  dis¬ 
tribution,  while  the  infrared  extinction  should  be  more 
strongly  affected  by  the  larger  particles  in  the  micron 
region  of  the  particle  size  distribution  [as  may  be  inferred 
from  figures  in  Nilsson  (1979)].  The  terms  accumulation 
and  coarse  modes  [Whitby  (1978)]  refer  to  two  peaks 
often  observed  at  two  size  ranges  in  the  particle  size  distri¬ 
bution.  Since  these  size  ranges  are  partly  influenced  by 
independent  processes,  the  two  modes  may  not  always  be 
closely  related  (Fitch  and  Cress  (1981)  and  Whitby 
(1978)].  Thus  one  might  expe-t  some  variation  in  the 
visible  to  infVired  extinction  relationship.  Jennings  dal. 
(1978)  discuss  the  effect  of  refractive  index  variations  on 
atmospheric  extinction. 

One  difficulty  inherent  in  the  aerosol  extinction  cal¬ 
culations  is  the  difficulty  in  assigning  reasonable  particle 
si/c  distribute  ns.  The  distributions  respond  to  a  number 
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of  environmental  factors,  both  macroscopic  and  micros¬ 
copic.  For  example,  Yue  and  Deepak  (1980)  discuss  the 
effects  of  sedimentation  rates  on  the  drop  size  distribution 
in  fogs. 

In  spite  of  the  uncertainties  in  particle  size  distribu¬ 
tion  and  refractive  mdex,  models  based  on  particle  size 
distributions  can  be  very  useful  Perhaps  the  most  com¬ 
monly  used  model  of  this  sort  :s  the  LOWTRAN  model. 
LOWTRAN  essentially  uses  four  haze  models  and  two  fog 
models  for  the  aerosol  portion  of  the  program.  The  haze 
models  are  based  on  an  assignment  of  the  most  likely  par¬ 
ticle  size  distribution  and  refractive  index  distribution  for 
four  conditions,  urban,  rural,  maritime,  and  tropospheric 
In  the  LOWTRAN  3B  model,  the  relative  particle  size  dis¬ 
tribution  a.aj  refractive  index  distribution  are  held  con¬ 
stant  within  each  model,  and  only  the  number  density  is 
allowed  to  vary.  Thus,  for  a  given  model,  say  rural,  the 
ratio  of  extinction  coefficient  at  two  given  wavelengths  is  a 
constant 

In  the  LOWTRAN  5  model  [Knetzys  etui.  (i980)j 
the  effect  of  relative  numidity  on  relative  particle  size  dis¬ 
tribution  and  refractive  mdex  is  included.  Thus  the  ratio 
of  infrared  to  visible  aerosol  extinction  coefficient  is  no 
longer  a  constant,  but  is  <.  function  of  relative  humidity 
The  effect  of  relative  humidity  is  much  stronger  in  the 
maritime  model  than  ir.  the  urban  or  rural  models  The 
analysis  in  this  report  will  include  some  discussion  of  how 
the  measured  infrared  to  visible  ratios  compare  with  the 
ratios  inherent  in  the  LOWTRAN  models. 


3.  DISCUSSION  OF  THE  MEASUREMENTS 

3.1  The  OPAQUE  Program 

The  OPAQUE  program  was  an  extensive  measure¬ 
ment  program  undertaken  by  several  nations  connected 
through  NATO  Measurements  were  taken  at  hourly 
intervals  over  a  two  year  period  at  several  stations 
throughout  Europe  The  program  is  well  documented  in 
Fenn  (1978).  The  analysis  leading  up  to  this  report  was 
begun  using  a  small  set  of  data  [Johnson,  et  at.  (1979b)). 
and  then  expanded  to  include  larger  portions  of  the  data 
base  Portions  of  the  Netherlands  data  set  were  chosen 
for  analysis  because  the  Netherlands  data  appeared  to  be 
the  most  complete  and  self  consistent  set  available  at  the 
time  The  present  report  is  based  upon  9  months  of  Neth¬ 
erlands  preliminary  release  data.  (See  Appendix  B.) 

The  measurements  were  made  by  the  Physics  Group 
of  .he  Physics  Laboratory  of  the  National  Defense 
Research  Organization-TNO,  at  Ypenberg  Air  Base 
ilanssen  (1981)]  The  Netherlands  OPAQUE  station  is 
located  at  the  Royal  Air  Force  Base  Ypenburg,  near  The 
Hague  and  Rotterdam.  It  is  near  the  coast,  but  is 
influenced  by  the  urban  industrial  environment  [Fenn 
(1978)1.  Figure  3-1  illustrates  the  location  of  the  Nether¬ 
lands  OPAQUE  station. 

The  OPAQUE  program  included  measurements  of 
visible  extinction  coefficient,  horizontal  and  vertical  illumi¬ 
nance.  path  radiance,  infrared  transmittance,  optical  and 
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infrared  turbulence,  atmospheric  turbidity,  standard 
meteorological  quantities,  and  aerosol  particle  size  distri¬ 
bution.  Of  these  data,  the  visible  extinction,  infrared 
transmittance,  and  meteorological  data  were  par ticularlv 
appropriate  for  this  study.  The  atmospheric  turbidity 
measurements  were  not  included  in  this  analysis  because 
in  the  fog,  when  the  infrared  signals  became  most 
significant,  the  pyrheiiometer  data  were  offscaie.  The 
aerosol  particle  data  would  be  of  particular  interest  too, 
however  the  instrument  was  not  yet  installed  during  the 
periods  analyzed  here  Aerosol  data  were  recorded  in  later 
seasons,  however  they  were  not  available  to  us  at  the  time 
of  this  analysis 

Thus,  after  some  preliminary  study  of  all  the  param¬ 
eters,  the  analysis  effort  was  directed  primarily  to  ihose 
parameters  summarized  in  Table  3.1. 

The  last  column  in  Table  3.1  lists  the  number  of 
available  data  points  per  measurement  period.  Each 
hourly  measurement  period  lasted  4  minutes 


3.2  Instrument  Discussion— Visible 

Two  measurements  of  visible  extinction  were  avail¬ 
able,  the  data  from  the  Eltro  transnnssometer  (1000m 
path  length)  converted  to  extinction  coefficients,  and  the 
scattering  coefficient  from  the  AEG  point  visibility  meter. 
The  goal  of  the  program  was  to  achieve  a  measurement 
accuracy  of  ±20%  in  scattering  coefficient  and  ±2'1. 
transmittance  For  the  purposes  of  this  analysis,  these 
numbers  will  be  used  as  the  most  probable  measurement 
accuracy.  (See  Appendix  B.)  The  ±2%  transmittance 
corresponds  to  an  uncertainty  in  extinction  coefficient  of 
about  5%  midrange,  increasing  to  30%  or  more  near  the 
limits  of  the  reported  transmittances.  In  prac.ice  it  is 
found  that  the  transnnssometer  indicates  somewhat  clearer 
air  than  the  nephelometer  does  at  the  high  transmittance 
end  of  the  scale.  Otherwise,  the  offsets  between  the 


instruments  are  not  particularly  systematic.  On  some 
days,  the  ,wo  instruments  are  in  agreement  within  close 
tolerance  all  day.  On  other  days  there  is  a  significant 
nearly  constant  offset,  while  on  others  there  is  a  variable 
offset. 

In  view  of  these  uncertainties  in  the  visible  extinc¬ 
tion.  we  initially  chose  for  infrared  comparison  two  sets  of 
5  days  in  which  the  two  visible  instruments  were  con¬ 
sistently  in  reasonably  close  agreement.  The  outlines  of 
the  prediction  scheme  were  developed  based  on  this  10 
aay  sample  The  analysis  was  then  extended  to  three  3- 
month  seasons  In  the  9  month  sample  analysis,  the  com¬ 
parisons  with  infrared  extinction  were  made  using  both 
the  nephelometer  data  and  the  transmissometer  data.  The 
nephelometer  data  have  yielded  slightly  more  consi  ,ent 
resuits  Because  of  this  consistency,  the  final  resuuc  in 
this  report  have  been  based  on  calculations  using  the 
nephelometer  data. 

3.3  Instrument  Discussion— Infrared 

The  Barnes  infrared  transmissometer  operates  over 
a  500  meter  path,  and  is  intended  to  achieve  a  ±2%  meas¬ 
urement  accuracy.  It  was  recognized  quite  early  in  the 
OPAQUE  program  that  there  were  difficulties  with  the 
standard  calibration  procedure  (Shand  (1978)].  As  a 
consequence,  a  different  calibration  scheme  was  evolved 
[Shettle  and  Fenn  (1978)  and  Kohnle  (1979)].  Basically, 
the  data  are  searched  for  a  clear  day  with  low  aerosol  and 
water  vapor  content,  i.c.  a  Jay  with  meteorological  range 
>20  km,  relative  humidity  ^80%.  dewpoint  temperature 
*$10CC,  and  no  piecipitalion.  (These  criteria  could  not 
always  be  met)  On  such  a  day.  the  aerosol  extinction 
should  ue  a  minor  portion  of  the  total  extinction,  and  the 
LOWTRAN  calculation  of  total  extinction  should  have 
optimal  accuracy  (Shelde.  1980)  Thus  a  calibration  con¬ 
stant  can  be  determined  by  comparing  the  measured 
transmittar.ee  on  this  clear  day  with  the  transmittance 


Table  3.1.  Measured  Variables  Used  for  Analysis 
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predicted  r  LOWTRAN  (including  the  aerosol  portion). 
The  uncertainty  in  the  calibration  constant  (determined  on 
a  clear  day)  results  in  an  uncertainty  in  the  calibration 
aerosol  extinction  from  other  days  which  becomes  quite 
small  in  haze,  mist,  or  fog.  This  uncertainty  is  discussed 
quantitatively  in  Section  5.2. 

4.  DATA  REDUCTION 

4.1  Derivation  of  Aerosol 
Extinction  Coefficients 

It  is  convenient  to  convert  the  visible  and  infrared 
data  to  aerosol  extinction  coefficients.  This  conversion 
expedites  both  the.  evaluation  of  the  variations  in  the 
measured  data  and  the  comparison  to  models. 

The  visible  data  were  supplied  in  the  form  of  total 
extinction  coefficients.  The  reported  extinction  coefficients 
range  from  .1  to  about  30  km'1,  while  the  molecular  or 
Rayleigh  component  is  about  .01  km  !  Since  this  Ray¬ 
leigh  component  is  much  le>s  than  die  instrumental  uncer¬ 
tainty,  the  measuieo  data  were  reated  as  aerosol  extinc¬ 
tion  coefficients.  Within  each  4-minute  measurement 
period,  a  beginning,  ending  maximum,  and  minimum 
value  were  given  All  four  measurements  were  used  in 
the  analysis,  but  the  results  presented  here  are  for  the 
average  of  the  beginning  and  ending  points. 

The  infrared  data  were  supplied  in  the  form  of  total 
transmittance  over  a  500m  path.  In  some  cases,  the  cali¬ 
bration  correction  mentioned  in  Section  >  had  already 
been  applied.  In  other  cases  the  correction  factor  was  sup¬ 
plied  and  the  data  were  corrected  accordingly.  The  aerosol 
transmittance  was  then  computed  using  the  following 
equation. 


Taer  -  ■  Tu„,(i)]  . 


In  this  computation,  the  water  vapor  (water  molecule)  and 
molecular  (air  molecule)  transmittances  (7//2o  and  T\fo,) 

were  computed  using  the  following  equations  from  Shettle 
(1978a). 


T,,/,.',)  =  1  -  ((•|+c2  /)exp((rfl+(/2-/)frf|  (4-2) 

T.i„,U)  =  <'+/■/  .  (4-3) 


These  equations,  derived  by  Shettle  on  a  partly  empirical 
and  partly  theoretical  basis,  represent  average  broadband 
transmittance,  derived  from  the  transmittances  of 
Lowtran.  as  in  Eq  (2-2)  The  appropriate  filter  responses 
of  the  instrument  were  used  in  Shettle's  derivation.  The 
accu'acy  of  Eqs.  (4-2)  and  (4-3)  relative  to  LOWTRAN  is 
better  than  1%  The  constants  of  the  equations  for  the 
Netherlands  instrument  [from  Shettle  (1978a)]  are  given 
in  Table  4.1. 

Since  Eq  (4-2)  includes  dewpoint  temperature  as  an 
input,  it  was  necessary  to  convert  the  reported  relative 


fable  4.1.  Constants  for  Molecular  and  Water  Vapor 
Transmittance  Computation  -  Netherlands. 
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humidities  to  dewpoint  temperatures.  This  was  done 
using  the  following  equation  froi.i  Shettle  (1978b). 


(19.772  i)  +  [ln(R///100)-(269.9+r)l  ,, 

d  =  19.772  -lnUW/100) 


Thus  the  results  of  Eq.  (4-4)  are  used  as  input  to  Eq.  (4- 
2).  Equation  (4-2)  is  used  along  with  Eq.  (4-3)  as  input 
to  Eq.  (4-1).  using  the  measured  temperatures  and  relative 
humidities.  Equation  (4-1)  yields  the  aerosol  transmit¬ 
tance,  which  is  then  converted  to  aerosol  extinction  using 
the  following  equation 


a„„(km  ') 


.5  (km) 


Thus  the  aerosol  extinction  is  derived  from  the  measured 
total  transmittance  and  the  computed  molecular  and  water 
vauor  transmittance  at  the  time  ot  measurement.  Since 
the  aerosol  extinction  is  reasonably  spectrally  invariant 
within  these  bands,  this  extinction  will  not  be  significantly 
range  dependent. 

For  the  band,  a  beginning  and  ending  meas¬ 

urement  were  supplied  for  each  4-minute  period.  Both 
were  used  in  the  analysis,  but  most  of  the  plots  utilize  the 
beginning  measurement. 

Another  quantity  used  for  comparison  to  the  aerosol 
extinction  was  the  seasonal  median  of  water  vapor  and 
molecular  extinction  effective  over  a  500m  path.  For  this 
quantity,  the  molecular  and  water  vapor  transmittances 
over  500m  were  computed  using  Eqs.  (4-2)  and  (4-3) 
(which  derive  from  LOWTRAN).  These  transmittances 
were  then  converted  to  effective  500m  extinctions  using 
Eq.  (2-5),  since  the  application  here  is  transmitted  radi¬ 
ance.  These  extinctions  showed  litt'e  variation  within  each 
season,  so  the  median  molecular  and  water  vapor  extinc¬ 
tion  was  used  in  various  figures  and  tables.  These  median 
500m  effective  extinctions  were  not  used  in  the  computa¬ 
tion  of  aerosol  extinction,  they  were  used  only  for  general 
comparison  to  the  aerosol  extinction. 


5.  DATA  ANALYSIS -DEVELOPMENT 
OF  THE  BIN  CONCEPT 

5.1  General  Parameter  Relationships 

The  general  relationships  between  the  infrared  aero¬ 
sol  extinction  and  the  other  parameters  may  be  con¬ 
veniently  illustrated  using  scatter  plots.  The  plots  in  Sec¬ 
tions  5.1  through  5.4  illustrate  the  data  from  one  season: 
June,  July  and  August  19'77.  Similar  plots  have  been  gen- 
era'ed  for  two  additional  season:,,  but  are  not  included  in 
the  report.  They  look  nearly  identical  to  the  plots  from 
Summer  ’77,  and  the  general  comments  in  this  section 
apply  equally  to  all  three  seasons. 

For  this  analysis,  the  data  during  rain  are  deleted. 
Rain  data  are  expected  to  be  characterized  by  less 
difference  between  infrared  extinction  and  visible  extinc¬ 
tion.  Since  the  scattering  properties  are  expected  to  be 
significantly  different  for  rain,  it  was  decided  to  omit  the 
rain  cases.  The  following  data  points  were  omitted:  any 
data  points  which  indicated  rain  during  the  measurement 
period  or  rain  during  that  hour;  and  any  points  for  which 
either  rain  measurement  was  not  recorded,  which  were 
thus  indeterminate. 

Figure  5-1  (a&b)  shows  scatter  plots  of  the  IR  aero¬ 
sol  extinction  coefficient  ( aa„ )  vs  the  visual  aerosol 
extinction  coefficient  (S),  in  the  3-5  and  8-12 (im  regions. 
The  median  effective  molecular  and  water  vapor  extinction 
(labeled  aSi0i+Hl o)  is  included  for  comparison.  In  Fig.  5-1 
there  is  a  fair  amount  of  statistical  scatter,  but  in  general 
the  IR  extinction  increases  as  the  visual  does  (for  this 
analysis,  since  we  are  dealing  exclusively  with  the  aerosol 
extinction  in  the  IR  bands,  it  is  convenient  to  refer  to  it 


simply  as  extinction).  The  squared  correlation  coefficient, 
r2,  between  log  aa„  and  log  5  (with  aat,  <  .01  omitted)  is 
26%  at  3-5 tim  and  21%  at  8-12 y.m.  Thus  the  relationship 
is  not  tight  enough  for  5  alone  to  be  an  accurate  predictor, 
yet  there  is  definitely  some  information  in  the  value  of  S. 
Much  of  the  variation  in  Fig.  5-1  is  expected  because  cf 
instrumental  uncertainty.  This  will  be  discussed  later  in 
this  section. 

Figure  5-2  (a&b)  shows  aerosol  extinction  ( aa„ )  at 
3-5 fim  vs  relative  humidity(R//)  and  aerosol  extinction 
(aa„)  at  8-12 nm  vs  RH.  The  !R  data  in  Fig.  5-2  have  a 
systematic  and  reasonable  behavior  with  respect  to  the 
relative  humidity.  Over  most  of  the  relative  humidity 
range,  the  IR  extinctions  are  low,  t.e.  mostly  between  01 
and  .1  km"1,  and  there  is  little  apparent  variation  with 
relative  humidity.  Above  humidities  of  about  94%  there 
are  a  large  number  of  high  infrared  extinctions  between  1 
and  10  km"1  or  higher.  High  infrared  extinctions  are 
expected  to  occur  in  the  presence  of  large  droplets  which 
should  normally  be  associated  with  nearly  saturated  air. 

It  is  interesting  to  note  in  these  figures  that  the 
highest  infrared  extinctions  were  associated  with  measured 
relative  humidities  near  95%.  These  data  points  often 
occurred  during  periods  when  the  measured  RH  values 
were  constant  for  sevetal  hours.  The  constant  values  vary 
from  one  day  to  the  next,  and  are  typically  values  such  as 
94.8%  and  95.3%.  The  constant  RH  signal  and  the  high 
values  of  IR  extinction  lead  us  to  believe  that  the  hygrom¬ 
eter  was  saturated,  and  that  the  true  relative  humidities 
were  probably  close  to  100%.  The  expected  instrumental 
uncertainty  is  about  5%  for  the  hygrometer.  (See  humi¬ 
dity  correction  in  Appendix  B) 


VISIBLE  SCATTERING  COEFFICIENT  (km1) 


VISIBLE  SCATTERING  COEFFICIENT  (km  ') 


Fig.  5-1.  Infrared  aerosol  extinction  vs  visible  scattering  coefficient 
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Fig.  5-2.  Infrared  aerosol  extinction  vs  relative  humidity. 


The  early  10-day  studies  yielded  no  systematic  rela¬ 
tion  between  infrared  aerosol  extinction  coefficient  and 
absolute  humidity.  Consequently  this  plot  was  not  gen¬ 
erated  for  the  3-month  intervals.  This  result  is  corro¬ 
borated  by  Bakker  (1981)  who  reported  an  analysis  of  a 
3-month  Fall  1978  sample  of  Netherlands  data.  He  notes 
"A  possible  interpretation  of  our  results  is  that  for  relative 
humidities  ^0.90  it  is  the  amount  of  liquid  water  [rather 
than  water  vapor  or  absolute  humidity],  present  in  the 
form  of  drops  in  rain  or  aerosois  in  fog  that  determines 
the  infrared  transmission  of  'he  atmosphere". 

Figure  5-3  illustrates  the  scatter  plot  for  the  visible 
extinction  coefficient  (5)  vs  the  relative  humidity.  The 
high  visible  extinctions,  like  the  high  IR  extinctions,  occur 
only  at  high  relative  humidities.  An  analysis  of  these  high 
extinctions  reveals  that  even  when  both  the  measured 
visual  extinction  and  relative  humidity  are  quite  high  (sug¬ 
gesting  the  presence  of  fog  and  large  droplets),  the  meas¬ 
ured  infrared  extinction  is  high  for  only  a  small  portion  of 
the  cases. 

Thus,  both  S  and  RH  have  predictive  value.  That 
is,  the  infrared  aerosol  extinctions  tend  to  increase  with 
the  visual  extinction  S.  The  highest  infiared  aerosol 
extinctions  are  found  at  high  relative  humidities,  but  even 
when  both  S  and  RH  are  high,  the  infrared  extinction  is 
not  always  high.  Thus,  the  goal  is  to  devise  the  best 
scheme  which  makes  use  of  the  information  in  both 
parameters.  Before  this  scheme  can  be  presented,  it  is 
important  to  evaluate  how  much  of  the  observed  variation 
in  the  parameter  relationships  might  be  attributed  to 
uncertainty  in  the  aerosol  extinction. 


5.2  Infrared  Aerosol  Extinction  Uncertainties 

In  evaluating  the  variation  in  the  infrared  aerosol 
extinction,  there  are  four  main  sources  of  uncertainty  to 
consider.  These  are  the  measurement  uncertainty,  the 
uncertainty  due  to  the  non-simultaneity  of  the  data,  the 


RELATIVE  HUMIDITY  (Percent) 


Fig.  5-3.  Visible  scattering  coefficient  (S)  vs  relative  humidity 
{RH),  summer  1977 
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uncertainty  associated  with  the  calibration  procedure,  and 
the  uncertainty  resulting  from  the  molecular  extinction. 
These  four  issues  wiil  be  b'ec.:ssed  in  the  following  para¬ 
graphs.  (See  Appendix  C ) 

The  projected  level  of  accuracy  for  the  infrared 
measurements  was  about  ±2%  in  total  transmittance. 
This  corresponds  to  a  variable  uncertainty  in  aerosol 
extinction,  as  shown  in  TaJe  5.1. 

Table  5.1.  Aerosol  Extinction  Uncertainty  for 
±2%  Uncertainty  in  Total  Transmittance. 

j  Aerosol  Extinction  Resulting  from  2#i  Change  .n  Toial  Transrpnl jneb  ' 
Aerosol  I  }-}uir  8-12am  , 
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1  were  computed 

using  the 

equation  (see  Appendix  C) 

m=d 

~  S'aacr  *  aMol+HJ 1  .1 

aaer 

e 

.5 

'  ±  U2) 

~am 

+„2o  (5-1) 

because  it  is  ±2%  of  total  transmittance.  For  this  compu¬ 
tation  the  median  molecular  plus  water  vapor  extinction 
for  Summer  ’77  was  used,  i.e.  .50  km-1  at  3-5 fim  and  .25 
km~!at  8-12 hm.  Note  that  the  uncertainty  corresponding 
to  ±2%  total  transmittance  is  primarily  a  function  of  the 
magnitude  of  the  aerosol  extinction,  as  may  be  seen  in 
Table  5.1.  The  uncertainty  depends  only  slightly  on  the 
magnitude  of  the  molecular  and  water  vapor  extinction, 
(or  on  the  relative  contribution  of  aerosol  to  total  extinc¬ 
tion)  resulting  in  slightly  different  values  in  Table  5.1  for 
the  3-5 band  and  the  8-12/a m  band.  Figure  5-4(a&b) 
shows  these  error  bars  superimposed  on  the  plots  of  aaer 
vs  S  (from  Fig.  5-1).  These  figures  also  show  the  pro¬ 
jected  ±20%  uncertainty  in  visible  scattering  coefficient, 
which  is  constant  in  S.  (See  Appendix  B.) 

In  Fig.  5-4,  note  that  most  of  the  variation  at  the 
clear  end,  i.e.  for  S  less  than  about  1  km'1,  can  be 
explained  simply  by  the  projected  infrared  measurement 
uncertainty.  For  conditions  where  5  is  greater  than  about 
1  km-1,  the  variation  in  the  infrared  extinction  becomes 
greater  than  the  projected  measurement  uncertainty.  Thus 
in  this  upper  region,  for  S>1  km"1,  there  is  a  significant 
amount  of  variation  which  might  be  explained  by  some 
predictor  such  as  S  or  RH. 

Although  the  plot  includes  only  infrared  aerosol 
extinctions  aar,  greater  than  .01  km"1,  the  data  set  includes 
a  significant  number  of  points  which  are  closer  to  0  or 
negative,  as  would  lie  expected  from  the  error  bars  illus¬ 
trated  in  Fig.  5-4.  Using  analysis  of  variance,  in  which  the 
range  in  S  was  divided  into  about  10  bins,  the  variance  in 
<xaet  within  each  of  these  bins  was  cumputr^.  For  this 


Fig.  5-4.  Uncertainly  in  infrared  aerosol  extinction  assoualed  with  a  2  uncertainly  in  measured  IR  transmittance  lcertic.il  bars) 
Horizontal  bar  represents  a  ±20^,,  uncctainty  in  visible  scattering  coefficient 
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analysis,  all  the  non-rain  data  were  used,  including  the 
negative  points.  The  variance  in  each  bin  was  then  com¬ 
pared  with  the  projected  uncertainty  at  the  midpoint  for 
each  bin. 

The  measured  variation  was  found  to  be  less  than  or 
equal  to  the  projected  measurement  uncertainty  for  all  the 
bins  with  S  less  than  about  1.  Here  the  variations  in  the 
measurements  of  a0fr  are  dominated  by  measuremen 
uncertainty,  and  should  not  be  predictable  from  other  con¬ 
siderations.  When  S'  becomes  greater  than  1  km-1,  how¬ 
ever,  the  observed  variation  is  much  more  than  the  meas¬ 
urement  uncertainty,  and  it  should  be  possible  to  explain 
or  predict  some  of  the  variation  in  aa„. 

The  second  source  of  uncertainty  mentioned  earlier 
is  the  uncertainty  due  to  the  non-simultaneity  of  the  data. 
At  the  Netherlands  station,  each  hour’s  data  taking 
occurrec  during  a  four  minute  period.  Since  the  1R  and 
visible  data  are  not  simultaneous,  this  could  contribute  to 
the  scatter  observed  in  Figs.  5-1  and  5-2.  Figure  5-5  com¬ 
pares  the  measurements  of  3-5 nm  aerosol  extinction 
which  were  taken  at  the  beginning  and  end  of  each  four 
minute  period.  Comparing  Fig.  5-5  (final  aL.,  vs  begin¬ 
ning  a <*,,)  with  Fig.  5-1  (a)  (beginning  aae/  vs  S)  shows 
that  the  variation  in  infrared  extinction  during  the  four 
minutes  is  much  less  than  the  observed  variation  in 
infrared  extinction  at  any  given  value  of  S.  The  squared 
correlation  constants,  r 2  are  26%  for  Fig.  5-1  (a)  and  74% 
for  Fig.  5-5. 

Thus  while  Fig.  5-4  illustrates  that  for  S  greater  than 
about  1  knr',  there  is  more  variation  than  can  be 
explained  by  measurement  uncertainty.  Fig.  5-5  illustrates 


Fip.  5-5.  final  beginning  infra-red  aerosol  exiinaion 
coefficient  al  3-5/im  fur  each  4  mmule  measuremenl  period 


that  this  variation  r  not  due  pnmarily  to  the  non¬ 
simultaneity  of  the  data. 

The  variations  in  S  during  the  four  minute  period 
were  also  studied  by  comparing  the  maximum  with  the 
minimum  S  and  the  beginning  with  the  final  S.  These 
variations  are  quite  small  relative  to  the  variation  in  Fig. 
5-1. 

The  third  source  of  uncertainty  is  the  calibration 
correction  scheme  mentioned  in  Section  3.3  All  of  the 
reported  data  have  been  corrected  by  a  calibration  constant 
determined  during  data  periods  with  high  visibility  and  low 
water  vapor  content.  During  such  periods,  the  aerosol 
extinction  should  be  low  so  the  LOWTRAN  estimated 
extinction  should  be  quite  accurate.  The  calibration  con¬ 
stant  is  determined  by  comparing  the  measured  value  with 
the  LOWTRAN  value  during  the  calibration  periods. 

If  there  were  no  1R  extinction  measurement  noise 
during  the  calibration  periods,  this  technique  would  force 
the  mean  IR  extinction  to  match  LOWTRAN  whenever  it 
is  clear.  In  actual  practice,  there  is  measurement  noise 
during  the  calibration  periods,  so  that  there  will  be  an 
uncertainty  in  the  calibration  constant.  Also,  there  will  be 
some  uncertainty  in  the  calibration  constant  due  to  uncer¬ 
tainty  in  the  visual  extinction,  temperature,  and  dewpoint 
used  to  compute  the  LOWTRAN  extinction  during  tne 
calibration  period.  Both  of  these  uncertainties  are  minim¬ 
ized  by  averaging  as  many  points  as  possible  to  determine 
the  calibration  constant.  In  general  it  is  safe  to  say  that 
when  the  measured  extinctions  are  high,  they  will  not  be 
significantly  affected  by  the  calibration  correction,  but  as 
conditions  approach  the  conditions  of  the  calibration 
period,  the  average  aerosol  extinctions  will  be  increasingly 
affected  by  the  calibration  correction. 

The  fourth  source  of  uncertainty  results  from  the 
subtraction  of  the  molecular  and  water  vapor  extinctions 
from  the  total  extinction  to  yield  the  aerosol  extinction. 
The  molecular  and  water  vapor  extinctions  are  functions 
of  the  temperature  and  the  relative  humidity.  An  uncer¬ 
tainty  of  1°  in  temperature  results  in  an  uncertainty  of 
.007  km"1  and  .014  km-1  in  the  final  infrared  aerosol 
extinctions  at  3-5 fim  and  8-12^m  respectively.  An  uncer¬ 
tainty  of  10%  in  relative  humidity  yields  an  uncertainty  of 
.018  km-1  and  .03  km'1  in  the  final  extinction  at  3-5 g/n 
and  8-12 nm  respectively.  These  uncertainties  will  have 
only  a  small  effect  on  the  low  extinctions  and  will  not 
significantly  affect  the  high  extinctions. 

The  net  effect  of  the  various  sources  of  uncertainty 
should  be  as  follows.  When  the  visual  extinction  is  low, 
i.e.  5<1  km1,  the  infrared  aerosol  extinction  is  generally 
low  (as  shown  in  Fig.  5-1).  In  this  reg..ne  the  measure¬ 
ment  uncertainty  should  yield  large  standard  deviations. 
The  use  of  the  calibration  constant  should  force  the  aver¬ 
age  infrared  aerosol  extinction  to  be  close  tc  LOWTRAN; 
how  close  the  average  is  will  depend  on  the  number  of 
measurements  averaged  to  determine  the  calibration  con¬ 
stant.  At  higher  values  of  infrared  extinction,  '  u-se 
uncertainties  should  have  little  effect,  but  theie  siould 
still  be  a  finite  but  small  uncertainty  due  to  ifie  i.?n- 
simultaneity  of  the  data.  The  effects  of  the  temperature 
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and  relative  humidity  uncertainties  on  the  infrared  extinc¬ 
tions  should  be  small  relative  to  these  effects. 

In  summary,  the  error  analysis  showed  that  when 
the  infrared  extinction  was  low,  most  of  the  variation 
could  be  attributed  to  measurement  uncertainty.  When 
the  infrared  extinction  was  high,  the  variation  could  not 
be  attributed  to  measurement  uncertainty  It  is  therefore 
reasonable  to  try  to  relate  the  variations  in  the  high 
infrared  extinctions  to  predictor  variations. 


5.3  Method  to  Sort  Infrared  Extinctions 

Taxing  into  consideration  the  conclusions  resulting 
from  the  scatter  plots,  summarized  in  Section  5.1,  and  the 
results  of  the  error  analysis,  summarized  in  Section  5.2, 
the  goal  is  to  devise  the  most  reasonable  means  of 
estimating  infrared  aerosol  extinction. 

For  these  data,  the  aerosol  extinction  is  usually 
much  less  than  the  median  molecular  and  extinctions.  (In 
Figs.  5-1  and  5-2  the  median  molecular  pius  water  vapor 
extinction  is  given.  The  molecular  and  water  vapor 
extinction  varies  with  temperature  and  humidity,  however 
within  each  season  its  standard  deviation  was  low.)  Con¬ 
veniently,  the  regime  which  is  dominated  by  measurement 
uncertainty  is  the  region  of  low  infrared  aerosol  extinc¬ 
tions,  in  which  the  uncertainty  in  aerosol  extinction 
represents  a  small  uncertainly  in  the  total  extinction.  That 
is.  when  the  aerosol  extinction  :s  low  the  variations  are 
not  meaningful  due  to  measurement  uncertainty  and  they 
are  not  significant  relative  to  the  molecular  and  water 
vapor  extinctions.  However,  when  the  aerosol  extinction 
is  high  the  variations  are  both  meaningful  (relative  to 
measurement  uncertainty)  and  significant  (relative  to  the 
molecular  and  water  vapor  extinctions). 

!n  view  of  these  considerations,  it  was  most  lmpor- 
"nt  to  try  to  develop  some  mechanism  for  sorting  out  the 
incidence  of  high  infrared  extinctions  (or  the  incidence  of 
high  infrared  to  visible  extinction  ratios).  If  the  low 
infrared  extinctions  could  be  isolated  by  some  criteria  in  a 
iower  bin.  they  could  be  treated  with  a  stochastic  model. 
The  remaining  data  would  go  into  an  upper  bin.  The 
upper  bin  would  hopefully  include  all  of  the  high  infrared 
extinctions.  These  extinctions  would  be  high  enough  that 
their  variations  become  more  significant  relative  to  the 
instrumental  noise,  and  they  would  also  be  high  enough 
relative  to  the  molecular  and  water  vapor  extinctions  that 
it  would  be  more  important  to  study  their  variations. 

The  really  high  infrared  aerosol  extinctions  are 
expected  to  occur  when  near  condensation  conditions 
result  in  suspended  water  droplets  of  size  about  equal  to 
the  wavelength  or  larger.  This  condition  is  defined  as  mist 
or  fog  depending  on  whether  the  visibility  is  greater  than  1 
km,  or  less  than  or  ecual  to  1  km  [McIntosh  (1963)  and 
Proulx  (1971)].  In  particular,  McIntosh  defines  mist  as: 
"A  state  of  atmospheric  obscurity  produced  by  suspended 
microscopic  water  droplets  or  wet  hygroscopic  particles 
The  term  is  used  for  synoptic  purposes  when  there  is  such 
obscurity  and  the  associated  visibility  exceeds  one  km:  the 
corresponding  relative  humidity  is  greater  than  about  95 


percent.  The  particles  contained  in  mist  have  diameters 
mainly  of  the  order  of  a  few  tens  of  microns.''  Thus  a  rea¬ 
sonable  approach  is  to  define  an  upper  bin  which  should 
include  ail  of  the  mist  and  fog  incidents 

It  should  oe  noted  that  in  the  U  S.,  the  term  mist  is 
popularly  usee  to  describe  drizz  e,  we  are  using  the  term 
not  to  mean  drizzle,  but  as  defined  above.  Also  it  should 
be  noted  that  it  is  Simmon  field  practice  to  report  this 
mist  condition  as  thin  fog,  judging  by  the  comments  of 
local  weather  offices  and  sample  v/eather  reports  from  the 
U.S.  and  F.urope  [see  Table  6-1  of  Duntley  eial.  (1973) 
and  Table  6-3  of  Johnson  et  al.  (1979a)]. 

Tnus  in  order  to  isolate  the  high  infrared  extinctions 
in  an  upper  bin.  a  reasonable  approach  is  to  attempt  to 
isolate  those  conditions  with  !ur0e  suspended  droplets,  n 
mist  and  fog.  For  this  data  base,  it  was  most  effective  to 
use  two  predictors  to  define  the  upper  bin  and  the  lower 
bin  First,  the  visible  extinction  is  required  to  be  greater 
than  1  km  1  (visibility  (  =  3  S)  <  3  km)  for  the  data  to  be 
included  in  the  upper  bin  (mist  and  fog  bin).  The 
S>1  km  'cutoff  was  chosen  because  it  is  a  reasonable 
threshold  for  mist,  this  is  where  the  .nstrumental  noise 
becomes  less  significant,  and  this  is  .he  threshold  above 
which  high  IR  aerosol  extinctions  begin  to  occur.  The 
additional  requirement  for  the  high  bin  data  ,s  that  RH 
must  be  greater  than  or  equal  to  about  94"  'The  cutoff 
used  was  log(l -RH/ 100)  =  -1  2.  or  RH  -9’  “%).  This 
threshold  limits  the  bin  to  near  condensation  cases  The 
choice  of  94%  as  opposed  to  95,,:  is  required  due  to  the 
uncertainty  in  the  hygrometer.  This  cutoff  allows  inclu¬ 
sion  of  the  very  high  infrared  extinctions  observed  at  rela¬ 
tive  humidities  above  94%  (Fig.  5-2)  Those  points  not 
included  in  the  upper  bin.  i.e.  those  with  S<1  km*'  or 
RH<  94%.  are  included  in  the  lower  bin.  (See 
Appendix  B.) 
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Figure  5-6  illustrates  the  separation  of  data  into  the  5-7 (a&b)  for  the  3-5^m  waveband,  and  Fig.  5-7 (c&d)  for 

lower  bin,  consisting  primarily  of  clear  to  haze  conditions,  the  8-12 /xm  waveband.  These  figures  are  histograms  of 

and  upper  bin,  consisting  primarily  of  mist  and  fog.  aaer.  They  show  the  number  of  occurrences  for  each 

Theoretically,  one  would  expect  this  sorting  mechanism  to  .05  km" 1  interval  in  aa„.  Figure  5-7 (a&c)  shows  the  aero- 

separate  the  bulk  of  the  low  infrared  extinctions  from  the  sol  extinctions  from  the  upper  bin,  defined  by  S>  1  km"1 

high  infrared  extinctions.  In  the  next  section  we  will  show  and\o%{\-RHHM)  ^-1.2.  The  remaining  aerosol  extinc- 

how  well  the  sorting  scheme  worked  tions  are  included  in  the  lower  bin,  shown  in  Fig. 

5-7(b&d),  and  defined  by  5<1  km-1  or  log(f- 
5,4  Results  of  Sorting  the  Data  /?/// 100)  >-1.2.  The  negative  values  shown  on  these 

plots  are  the  result  ol  the  measurement  uncertainties  dis- 
The  results  of  sorting  the  data  are  shown  in  Fig.  cussed  in  Section  5.2. 
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Fig.  5-7.  Hisiograms  of  infrared  aerosol  extinction,  3-5/jm  upper  bin  (Mist!.  3-5jum  lower  bin  (Haze).  S-i2/um  upper  bin. 

and  8-!2/um  lower  bin. 


From  these  plots,  one  can  see  that  the  method  was 
reasonably  successful  in  sorting  out  the  high  infrared 
extinctions.  Nearly  all  the  high  infrareo  extinctions  are 
included  in  the  upper  bin.  The  bins  do  not  appear  to  be 
improved  by  slightly  changing  the  thresholds.  If  the 
cutoffs  were  raised  in  either  parameter,  (i.e.  if  one  used  a 
higher  5  or  higher  RH  cutoff),  then  too  many  high 
infrared  extinctions  were  included  m  the  lower  bin;  if  the 
cutoffs  were  lowered  in  either  parameter,  then  more  low 
extinctions  were  included  in  the  upper  bin.  Also,  we 
found  the  use  of  both  parameters  as  predictors  necessary. 
If  either  parameter  were  used  alone,  then  a  large  number 
of  low  extinctions  were  thrown  into  the  upper  bin. 

Note  that  the  upper  bin,  while  including  most  of  the 
high  extinctions,  also  includes  a  great  many  low  extinc¬ 
tions.  As  shown  in  later  sections,  the  resulting  large  vari¬ 
ation  in  auc,  in  the  upper  bin  is  not  due  to  instrumental 
uncertainty,  but  appears  to  be  the  result  of  variations  in 
the  aerosol  properties. 

5.S  Application  of  the  Sorting 
Mechanism  to  Other  Seasons 

Although  it  is  reasonable  to  expect  that  most  of  the 
high  IR  extinctions  would  be  associated  with  mist  and  fog 
episodes,  it  is  important  to  explore  the  effectiveness  of  the 
S  and  RH  cutoff  criteria  with  an  extended  data  base.  Two 
other  3-month  seasons  were  chosen.  June,  July  and 
August  of  1978;  and  December  1977,  January  1978,  and 
February  1978.  We  designated  the  three  3-month  seasons 
Summer  ’77,  Summer  ’78,  and  Winter  ’77. 

The  scatter  plots  of  the  other  two  seasons  looked 
very  much  like  the  scatter  plots  from  Summer  ’77  illus¬ 
trated  in  Section  5.1.  The  sorting  mechanism  described  in 
previous  sections  were  applied  to  the  other  two  seasons, 
and  worked  very  well.  In  each  season,  the  lower  bin 
includes  very  few  occurrences  of  high  infrared  extinction. 
The  data  for  all  three  seasons  have  been  converted  to 
cumulative  frequency  plots,  as  shown  in  Fig  5-8  (a-f). 
These  cumulative  frequency  plots  show,  for  each  value  of 
a0fr,  the  percentage  of  car  vith  infrared  aerosol  extinc¬ 
tion  less  than  or  equal  tu  that  value.  (Note  the  scale 
change  at  aa„«. 3.)  In  all  cases,  the  low  extinctions  are 
much  more  prevalent  in  the  lower  bin,  and  the  high 
extinctions  are  essentially  solated  in  the  upper  bin. 

The  median  values  of  infrared  aerosol  extinction 
and  the  standard  deviations  for  each  bin  are  summarized 
in  Tables  5.2  and  5.?.  The  median  values  correspond  to 
the  50%  values  in  Fig.  5-8.  In  these  tables  the  seasonal 
changes  are  much  less  than  the  differences  between  the 
upper  and  lower  bins.  In  all  cases  the  seasonal  changes 
are  .’•ell  within  the  standard  deviations  of  the  bins  These 
standard  deviations  are  strongly  influenced  by  a  few  high 
values,  however  they  give  a  general  idea  of  the  variability. 
In  the  lower  bin,  both  the  median  values  and  the  standard 
deviations  are  in  all  cases  very  low  relative  to  the  median 
molecular  and  water  vapor  extinctions. 

For  the  upper  bin,  the  median  values  are  in  all  cases 
(i.e.  all  seasons  and  filters)  much  higher  than  the  medians 


Table  5.2.  Bin  Statistics,  Infrared  Aerosol  Extinction 
( acer  at  3— Sum) 
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Table  5.3.  Bin  Statistics,  Infrared  Aerosol  Extinction 
(aaer  at  8-12 fj.m) 


Season 
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Aerosol  Extinction 
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03 

lor  the  lower  bin.  More  significantly,  the  standard  devia¬ 
tion;  are  quite  large  in  the  upper  bin,  even  relative  to  the 
magnitude  of  the  molecuiar  extinction.  Thus  in  the  upper 
bin  the  variations  become  more  significant,  and  worth 
looking  at  in  greater  detaii. 

Thus  by  using  two  predictors  for  son.ng,  the  visible 
extinction  and  the  relative  humidity,  one  can  create  two 
bins:  a  bin  which  consists  almost  entirely  of  low  values  in 
infrared  extinction,  and  a  smaller  high  bin  which  includes 
nearly  all  the  high  values  of  infrared  extinction  as  well  as 
a  significant  number  of  low  values.  This  sorting  mechan¬ 
ism  has  been  applied  to  two  summer  3-month  seasons  and 
one  winter  3-month  season,  and  seems  to  apply  well  to  all 
three  seasons.  This  is  a  useful  predictive  scheme,  in  that 
it  isolates  the  vast  majority  of  the  high  extinctions  in  one 
bin. 


6.  DATA  ANALYSIS-CHARACTERISTICS 
OF  THE  BINS 

6.1  Characteristics  of  the  Lower  Bin 

The  lower  bin  or  haze  category  is  defined  as  the  set 
of  points  with  the  visual  extinction  coefficients  less  than  or 
equal  to  1  (km-1),  or  with  log (1 — RH/100)  greater  than 
-1.2  (i.e.  relative  humidity  <94%).  That  is,  the  lower  bin 
includes  the  lower  end  of  the  visual  extinction  and  relative 
humidity  ranges.  This  bin  contains  the  vast  majority  of 
t  te  points;  92%  and  91%  in  the  summers  of  '77  and  ’78, 
and  72%  in  the  winter. 

The  infrared  extinctions  are  consistently  low  in  this 
bin.  The  observed  variation  is  commensurate  with  a 
measurement  uncertainty  of  ±2%  transmittance  as  shown 
in  Table  6.1.  The  observed  variation  in  column  3  (STD  in 
aa„)  is  the  standard  deviation  in  the  infrared  aerosol 
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CUMULATIVE  FREQUENCY  (%)  CUMULATIVE  FREQUENCY  (%)  CUMULATIVE  FREQUENCY  (%) 


Table  6.1.  Obsei  .'•id  Variation  in  Infrared  Aerosol 
Extinction  (a„r)  for  Lower  Bin  Compared 
with  Projected  Measurement  Uncertainty  (±2%T). 
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extinction  ow  in  the  lower  bin.  For  two  of  the  seasons, 
these  standard  deviations  are  strongly  influenced  by  1  or  2 
very  high  values  which  may  be  spurious.  Column  4  shows 
the  standard  deviations  with  these  1  or  2  exceptional 
points  deleted.  (For  Summer  ’77,  Summer  ’78,  and 
Winter  ’77,  the  number  of  points  deleted  were  respec¬ 
tively:  2  out  of  173.2,  0  out  of  1240,  and  1  out  of  1112.) 
In  column  5  the  2%  uncertainty  in  transmittance  has  been 
converted  to  an  uncertainty  in  aaer  at  the  median 

in  Table  6.1,  the  observed  variation  in  aerosol 
extinction  in  the  lower  bin  is  close  to  the  projected  meas¬ 
urement  uncertainty.  This  means  that  whatever  variation 
may  have  been  caused  by  changes  in  particle  size  distribu¬ 
tion  or  refractive  index  distribution  of  the  aerosol,  these 
variations  are  not  much  larger  than  the  projec'ed  measure¬ 
ment  uncertainty.  One  should  not  expect  to  be  able  to 
explain  much  of  the  real  variation  using  other  predictors 
because  the  real  variations  should  be  masked  by  measure¬ 
ment  uncertainty. 

Although  systematic  relations  between  the  aerosol 
extinctions  and  available  predictors  were  not  expected  (on 
the  basis  of  the  projected  infrar’d  measurement  uncer¬ 
tainty),  one  can  not  be  certain  ot  me  character  and  magni¬ 
tude  of  the  actual  infrared  measurement  uncertainty.  For 
this  reason,  some  study  was  made  of  the  parameter  rela¬ 
tionships  in  the  lower  bin.  The  ratio  of  infrared  aerosol  to 
visible  extinction  ( aarr/S )  was  evaluated  and  found  to  be 
more  variable  than  the  infrared  aerosol  extinctions.  The 
character  of  the  aaer/S  ratio  changes  appears  to  indicate 
that  much  of  the  ratio  variation  in  the  lower  bin  is  due  to 
measurement  uncertainty  in  the  visible  extinctions  at  the 
dear  end. 

Aiso  an  attempt  was  made  to  further  separate  the 
lower  bin  cases  into  three  bins  on  the  basis  of  relative 
humidity  and  visible  extinction.  This  method  did  not  yield 
useful  infrared  aerosol  extinction  separations.  This 
appears  to  support  the  expectation  that  the  real  variations 
in  the  lower  bin  are  most  probably  masked  by  measure¬ 
ment  uncertainty. 

If  further  precision  is  desired  in  this  lower  bu  it 
would  be  best  to  base  the  model  on  measurements  taken 
with  a  longer  path  length  instrument  in  the  infrared  as 
well  as  in  the  visible,  to  avoid  much  of  the  contamination 
caused  by  measurement  uncertainty.  For  most  purposes 
however,  it  is  not  necessary  to  resolve  the  variations  in 


the  lower  (haze)  bin,  since  the  aerosol  extinctions  are  con¬ 
sistently  low  relative  to  the  molecular  and  water  vapor 
extinctions. 

6.2  Estimation  of  Infrared  Extinction  for  Haze  Cases 

This  section  describes  how  to  use  the  statistics  in 
the  lower  bin  to  provide  an  estimation  of  infrared  aerosol 
extinction  for  at  least  this  location.  (The  statistics  may  be 
site  specific,  however  the  general  conclusions  may  apply 
more  extensively.)  As  described  in  the  previous  section,  <t 
was  found  that  the  infrared  aerosol  extinct'on  a0„,  rather 
than  the  extinction  ratio  aaer/S ,  was  better  behaved  in  the 
lower  bin.  As  can  be  seen  in  the  histograms  for  the  lower 
bin  from  Summer  ’77,  in  Fig.  5-7(b&d),  the  measured 
distribution  of  aaer  in  the  lowei  bin  is  not  a  normal  distri¬ 
bution.  Because  of  this,  cumulative  frequency  tables  such 
as  Table  6.2  are  used  to  describe  the  observed  variance  in 
infrared  extinction.  In  Table  6.2  the  two  summers  have 
been  combined,  to  yield  the  most  representative  summer 
distribution. 

This  table  may  be  interpreted  as  follows.  The  most 
typical  value  of  measured  infrared  aerosol  extinction  for 
this  location  is  given  by  the  50%  value.  The  extinction 
that  was  exceeded  only  5%  of  the  time  is  given  by  the  $5% 
value.  Since  much  of  this  measured  variation  may  be 
attributed  to  measurement  uncertainty,  one  cannot  say 
that  the  actual  variation  iri  the  extinction  is  as  large  -s 
implied  in  Table  6.2,  hov  ever  it  is  safe  to  say  that  the 
actual  variation  in  the  extinction  is  less  than  or  equal  to 
the  measured  variation.  The  negative  values  in  Table  6.2 
are  an  artifact  of  this  measurement  uncertainty 

Table  6.2  demonstrates  what  we  have  stated  before, 
that  in  non-fog  cases  the  aerosol  extinction  is  consistently 
low.  The  95%  value  is  only  about  1  km"1.  Even  the  98% 
values  are  in  most  cases  much  lower  than  the  molecular 
and  water  vapor  extinction.  As  long  as  the  visible  extinc¬ 
tion  is  less  than  1  or  the  relative  humidity  is  less  than 
94%.  aerosol  extinctions  which  approach  molecular  and 
water  vapor  extinction  are  very  rare  events. 

The  expected  error  in  estimation  for  this  location 
may  be  illustrated  by  comparing  the  measured  va'ues  for 
one  season  with  the  values  estimated  on  the  basis  of 
Table  6.2.  The  50%  values  from  Table  6.2  were  used 


Table  6.2.  Cumulative  Frequency,  Aerosol  Extinction 
(aoer)  Lower  Bin  (Haze) 
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along  with  the  measured  a values  to  compute  the 
cumulative  frequency  for  the  square  of  the  prediction 
error  (the  square  of  the  difference  between  and  the 
50%  value)  The  cumulative  frequency  for  the  error 
squared  is  shown  in  Fig  6-1  In  these  plots,  the  x  axis 
ranges  from  0  to  (u.w„,  +  a„2 0)2  (rounded  to  a  convenient 
number)  using  the  median  aWL  +a/,l0  from  Summer 
'll  This  scale  was  chosen  because  the  absolute  value  of 
the  error  is  of  less  importance  than  the  magnitude  of  the 
error  relative  to  the  molecular  and  water  vapor  extinction. 
The  median  error2  values  are  listed  on  the  plot.  These 
plots  illustrate  that  the  estimation  errors  are  quite  low 
For  both  niters,  the  estimation  errors  are  much  less  than 
the  molecular  extinction  at  least  99%  of  the  time  except 
for  cxaer  (8— 12/um)  in  the  winter,  when  the  percentage 
drops  to  97%  due  to  the  low  molecular  extinction. 

It  should  be  noted  that  the  aerosol  extinction  shows 
very  little  spectral  change  within  these  wavebands  (based 
on  LOWTRAN5  and  the  LOWTRAN  aerosol 
specifications).  Thus  the  %'alues  in  the  cumulative  fre¬ 
quency  tables  should  apply  reasonably  well  to  most 
wavelengths  within  the  two  wavebands.  To  obtain  total 
extinction,  some  technique  such  as  the  LOWTRAN  model 
must  be  used  to  calculate  the  molecuiar  and  water  vapor 
extinctions  for  the  appropriate  wavelength,  since  these  are 
highly  spectrally  dependent.  These  derived  model  molecu¬ 
lar  and  water  vapor  extinctions  must  be  added  to  the 
estimated  aerosol  extinction  to  obtain  the  estimated  total 
extinction  coefficient. 


6.3  Comparison  of  Measured  and 
Model  Haze  Aerosol  Extinctions 

As  an  integral  part  of  the  program  LOWTRAN5, 
the  clear  to  haze  aerosol  extinction  is  calculated  from  Mie 
theory  assuming  several  modelled  particle  size  distribution 
and  refractive  index  combinations,  to  represent  rural  con¬ 
ditions.  urban,  etc.  Within  each  model,  the  relative  size 
distribution  and  refractive  index  are  allowed  to  vary  with 
relative  humidity,  and  the  number  density  is  normalized 
by  the  visual  extinction.  As  a  result,  within  each  model, 
the  infrared  to  visible  ratio  is  a  1  unction  omy  of  relative 
humidity. 

In  another  model  by  Huschke  (1976),  the  8—12/am 
aerosol  extinction  is  derived  from  t  visible  extinction  and 
the  ratio  of  infrared  aerosol  extinction  to  visible  extinc¬ 
tion.  This  ratio  remains  constant  for  relative  humidities 
•ess  than  or  equal  to  50%,  and  the  ratio  increases  to  a  ratio 
of  1.0  at  a  relative  humidity  of  100%. 

Table  6.3  compares  the  measured  low  bin  extinc¬ 
tions  with  the  aerosol  extinctions  derived  from  the 
LOWTRAN5  and  Huschke  models.  In  this  table,  the 
measured  Maximum,  Median,  and  Minimum  values  are 
represented  by  the  98%,  50%,  and  2%  cumulative  fre¬ 
quency  measured  extinctions.  The  model  Maximum, 
Median,  and  Minimum  values  are  the  model  values  com¬ 
puted  from  the  98%  cumulative  frequency  value  of  S  and 
RH= 99%;  from  the  median  S  and  RH.  and  from  the 
minimum  S  and  RH.  The  values  in  Table  6.3  are  also 
plotted  in  Fig.  6-2. 
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Fig.  6-1.  Cumulatixe  Irsqucno  of  squared  error  ol  estimate  for  aerosol  extinction  at  3-5/um.  lower  bin  (Haze),  three  3-month  seasons 
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In  Table  6.3  and  Fig.  6-2,  the  median  aerosol  extinc¬ 
tion  is  generally  slightly  lower  than  predicted  by  the 
LOWTRAN  rural,  urban,  and  maritime  models,  but 
slightly  higher  than  predicted  by  the  LOWTRAN  tropos¬ 
pheric  and  Huschke  models.  The  models  all  fall  well 
within  the  range  of  measured  extinctions. 

6.4  Summary  for  I-ower  Bin  (Hazel  Cases 

It  was  found  that  if  either  the  visible  extinction  is 
less  than  or  equal  to  1  km'1  or  the  relative  humidity  is  less 
than  94%,  the  data  for  this  location  may  be  classed  in  a 
low  bin  containing  primarily  clear  to  haze  conditions.  In 
this  bin  the  infrared  aerosol  extinctions  compare  reason¬ 
ably  well  to  model  predictions.  The  aerosol  extinctions  are 
consistently  low  relative  to  molecular  and  water  vapor 
extinctions  (thus  they  are  a  small  portion  of  the  total 
extinction)  and  may  be  represented  by  cumulative  fre¬ 
quency  of  occurrence  tables. 


Table  6.3.  Comparison  of  Measured  and  Modelled  Aerosol 
Extinctions  for  Lower  Bin  (Haze). 
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6.5  Characteristics  of  the  Upper  Bin 

The  upper  bin  was  designed  to  include  all  the  mist 
and  fog  data.  It  is  defined  as  the  set  of  points  with  the 
visual  extinction  coefficient  S  greater  than  1  km  1  and 
\og(l-RH/\00)  less  than  or  equal  to  -1.2  (i.e.  relative 
humidity  ^94%).  That  is,  the  upper  bin  includes  the 
points  which  have  both  high  visual  extinction  and  high 


relative  humidity.  The  bin  contains  about  8%  of  the 
points  in  the  summer  and  about  22%  of  the  points  in  the 
winter.  It  contains  nearly  all  the  high  infrared  aerosol 
extinctions,  as  well  as  a  significant  number  of  low  infrared 
aerosol  extinctions. 
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Fig.  6-2.  Comparison  of  m:asured  and  model  infrared  aerosol  extinctions,  lower  bin  (Ha/e) 


A  test  of  the  effectiveness  of  the  bin  thresholds  may 
be  performed  by  evaluating  the  ratio  of  infrared  aerosol 
extinction  to  visible  xtinction,  <x„JS.  This  ratio  is 
expected  to  approach  values  near  1  only  in  the  presence  of 
large  droplets  (of  size  ~10 fj.m  or  more)  such  as  may 
occur  in  mist  or  fog.  The  histograms  of  aafr/S  for  the 
upper  bin  from  Summer  ’77,  shown  in  Fig.  6-3,  illustrate 
that  there  are  a  significant  number  of  cases  with  ratios 
near  1  in  this  bin.  By  contrast,  the  lower  bin  has  some¬ 
what  high  ratios  only  when  instrumental  noise  in  S 
becomes  large.  In  Fig.  6-3 (a),  there  are  several  values 
with  ratios  above  1,  however  this  is  consistent  with  the  3- 
5 fitn  LOWTRAN  estimate  for  Radiation  Fog  Model  2 
[Shettle  and  Fenn  (1979)],  which  is  1.45. 

An  interesting  feature  cf  these  histograms  is  the 
wide  range  in  the  magnitude  of  the  a„e,/S  ratio,  and  the 
prevalence  of  low  ratios.  The  high  ratios  are  distributed 
over  the  whole  range  in  5,  .n  the  upper  bin.  Calculations 
utilizing  the  ending  3-5 /xm  point  (T5),  and  the  beginning 
and  ending  visible  points  (SMAX  and  SMIN)  showed 
similar  distributions. 

Like  the  aae,/S  ratio,  the  aerosol  extinction  aaer 
varies  over  a  wide  range  in  the  upper  bin  (ref.  Fig.  5-8). 
The  distribution  in  aerosol  extinctions  is  summarized  in 
Table  6.4.  This  table  lists  the  percentage  of  cases  which 
exceed  three  thresholds:  <xacr  >  .1  km-1,  aae,  >  median 
<*moI+h2 o  for  each  season  as  listed  in  Tables  5.2  and  5.3, 
and  >  1.0  km'1. 

Note  that  even  in  this  mist  bin,  the  infrared  aerosol 
extinction  is  greater  than  the  molecular  and  water  extinc- 


Table  6.4.  Aerosol  Extinction  Values  in  Upper  Bin. 


Waveb«nd 

— 

Season 

>  1 

%  of  C*ses  with  <xQ(f 
>  Me1'  aMoI  +  H2 0 

>  1 

3*5/x  m 

Summer  ’77 

63 

24 

17 

Summer  ’78 

38 

3 

2 

Winter  ’77 

62 

28 

17 

8-12  fim 

Summer  ’77 

54 

23 

8 

Summer  ’78 

41 

3 

2 

Winter  '77 

43 

34 

13 

tions  only  about  20%  of  the  time,  and  greater  than 
1.0  km-1  only  about  10%  of  the  time,  on  the  average. 
These  high  extinctions  occur  often  enough  to  be  important 
in  many  applications,  however  one  certainly  cannot 
assume  that  the  infrared  aerosol  extinction  will  be  high  at 
all  times. 

Quite  a  bit  of  analysis  has  beer  done  to  determine 
whether  the  variations  in  infrared  aerosol  extinction  in  the 
upper  bin  could  be  related  to  some  predictor.  There  is  a 
tendency  for  the  infrared  aerosol  extinction  a ^  to 
increase  along  with  the  visible  extinctions.  As  a  result, 
the  aa„/S  ratio  is  somewhat  less  variable  than  the  extinc¬ 
tion  aa„.  There  are  a  number  of  statistical  tests  which 
illustrate  this,  the  most  important  being  that  the  final  error 
of  estimate  is  reduced  if  the  estimation  is  based  on  the 
median  a„„/S  ratio  rather  than  the  median  extinction 
ocaer. 

Even  though  the  ratio  aat,/S  is  better  behaved  than 


IR  AEROSOL  EXTINCTION  C->m)  /  IR  AEROSOL  EXTINCTION  (8-12/im)  / 

VISIBLE  SCATTERING  COEFFICIENT  VISIBLE  SCATTERING  COEFFICIENT 

Fig.  6-3.  Histograms  of  ihe  ratio  of  infrared  aerosol  extinction  lo  visible  scattering  coefficient,  upper  bin  (Mist)  summer  1977 
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the  extinction  aa„,  the  ratio  is  quite  variable  (as  noted 
earlier).  By  utilizing  the  ratio  aa„/S ,  we  are  essentially 
removing  the  effects  of  changes  in  the  magnitude  of  parti¬ 
cle  number  density  Changes  in  relative  particle  size  dis¬ 
tribution  or  refractive  index  would  be  expected  to  affect 
the  a  a,, /S  ratio. 

An  analysis  was  made  of  the  relationship  between 
the  a„(r/S  ratio  and  a  variety  of  parameters,  such  as  wind 
speed  and  wind  direction,  and  found  no  consistent  rela¬ 
tionships.  Particularly  the  relationship  between  the  aae,/S 
ratio  and  the  visible  extinction  S  was  studied,  however 
analysis  or  variance  showed  that  what  little  relationship 
existed  between  aae,/S  and  S  for  the  Summer  'll  data 
was  reversed  in  the  Summer  ’78  data.  Figure  6-4  shows 
aaer/S  vs  S  for  Summer  ’77  Some  sample  error  bars 
have  been  added  to  Fig.  6-4(a).  The  error  in  aae,/S  was 
computed  assuming  the  errors  in  aaer  and  in  5  are 
uncorrelated.  The  errors  for  Fig.  6-4(b)  arc  of  similar 
magnitude. 

In  Fig.  6-4,  the  high  uai.rfS  ratios  are  distributed 
over  the  whole  range  in  5  m  the  upper  bin  -  they  are  not 
associated  only  with  the  highlv  turbid  condition;.  The 
significance  of  this  is  that  no  matter  what  S  threshold  is 
used  to  define  fog,  the  e  will  still  be  a  mixture  of  high  and 
low  ratios.  Fog  is  internationally  defined  to  have  meteoro¬ 
logical  range  <1  km,  or  5  ^3.9  km'1.  Looking  at  the 
cases  with  S>4  km"'  in  Fig.  6-4,  one  can  see  that  at  least 
half  the  points  still  have  low  ratios.  Conversely,  in  the 
1<S<4  km"'  range,  which  is  defined  not  to  be  fog,  there 
are  still  a  significant  number  of  cases  with  high  ratios. 
That  is,  in  both  mist  and  fog,  the  infrared  to  visible  ratio 
is  often  high,  but  even  in  fog  it  is  also  often  low 


Some  of  the  variation  in  the  infrared  to  visible  ratio 
in  the  upper  bin  could  potentially  be  related  to  relative 
humidity.  We  do  not  have  a  precise  measure  of  relative 
humidity.  The  relative  humidity  measurement  accuracy  is 
about  ±6%,  judging  by  the  fact  that  there  were  a  number 
of  cases  with  constant  relative  humidity  of  about  94% 
associated  with  high  infrared  to  -'isible  ratios  (which  indi¬ 
cates  the  true  relative  humidity  was  most  probably  about 
100%)  Come  of  the  data  in  the  upper  bin  were  almost 
undoubtabiy  saturated,  but  others  may  not  have  been,  and 
we  have  no  way  to  separate  them.  Future  measurements 
using  more  precise  instrumentation  to  determine  relative 
humidity  might  allow  one  to  partially  separate  the  high 
ratios  from  the  low  ones. 

It  should  also  be  noted  that  some  of  the  variation  in 
the  aaf,/S  ratio  may  be  related  to  aerosol  type.  Some  of 
these  aerosols  may  have  been  advection  fogs,  and  others 
may  have  been  radiation  fogs.  The  advection  fogs  are 
normally  associated  with  large  droplets  and  high  infrared 
to  visible  ratios,  while  the  radiation  fogs  are  normally 
associated  with  small  droplets  and  lo  ver  infrared  to  visible 
ratios.  Examination  of  the  synoptic  and  mesoscale 
meteorology  would  be  useful  in  exploring  this  possibility. 

6.6  Study  of  the  IR  Aerosol  Extinction  Persistence 
in  the  Upper  Bin 

In  the  previous  section,  it  was  found  that  even  in 
the  mist  bin,  the  infrared  aerosol  extinction  remained 
lower  than  the  molecular  and  water  vapor  extinction  about 
80%  of  the  time.  The  incidence  of  the  high  extinctions 
did  not  appear  to  be  related  to  available  predictors.  It  is 


Fig.  6-4.  Ratio  of  infrared  aerosol  extinction  lo  visible  extinction.  aa„/S.  vs  visible  extinction  S'tkm  '). 

upper  bin  (Mist),  summer  1977. 


important  to  know  how  stable  the  infrared  conditions 
were,  i.e.  how  long  the  infrared  aerosol  extinctions  per¬ 
sisted. 

in  the  absence  of  significant  periodic  variations,  the 
correlation  coefficient  between  measurements  of  the  same 
parameter  at  times  0  and  /  would  be  expected  to  decrease 
exponentially  with  time  /  according  to  the  equation 

/-</)  =  /-(We  (6-1) 

where  a  is  the  relaxation  constant,  and  r( 0)  is  the  correla¬ 
tion  coefficient  for  Atime=0.  The  value  of  /TO)  is  related 
to  measurement  uncertainty.  (See  Appendix  B.) 

This  type  of  calculation  is  called  a  persistence  calcu¬ 
lation.  Note  that  these  persistence  calculations  are  not  cal¬ 
culations  of  the  correlation  between  the  visible  and 
infrared,  nor  are  they  frequency  of  occurrence  calcula¬ 
tions  Persistence  calculations  relate  to  the  question  of 
how  stable  a  given  parameter  is  with  the  passage  of  time. 
See  Brooks  and  Carruthers  (1953) 

For  this  study,  calculations  were  made  of  the  time 
variations  in  the  infrared  as  well  as  in  the  visible.  Two 
correlations  were  run  in  the  3-5 (im  waveband,  one  in  the 
8-12 ix m  waveband,  and  two  in  the  visible.  In  the 
3-Sfj.ni  band  and  the  visible,  the  measurements  taken  at 
the  beginning  of  each  measuremeni  period  were  correlated 
with  those  at  the  end  of  the  4-minute  period,  and  were 
correlated  with  the  values  taken  at  the  next  measurement 
period  60  minutes  later  The  two  resulting  correlation 
coefficients  may  be  represented  by  /T4)  and  a  (60).  For 
the  8-12 band,  only  r (60)  could  be  computed. 

The  squared  correlation  coefficients  r(4)  and  r(60' 
are  listed  in  Table  6.5.  These  statistics  are  based  on  In  « 
or  In  5.  The  r 2  values  for  60  minutes  are  fairly  low  in 
both  the  infrared  bands,  but  somewhat  higher  in  the  visi¬ 
ble  band 

The  persistence  correlations  for  the  3-5  ///» 
waveband  are  plotted  in  Fig.  6-5.  Note  that  the  intercepts 
are  slightly  less  than  I  as  expected.  The  values  of  /TO)  for 
the  3-5/x///  band  are  listed  in  column  2  of  Table  6  6. 
(These  values  are  approximate,  due  to  the  uncertainties  in 
the  extrapolation.)  The  r(0)  coefficients  may  be  converted 
to  estimates  of  standard  errors  using  the  following  equa¬ 
tion,  which  applies  to  normal  distributions. 


Table  6.5.  Squared  Persistence  Correlation  Coefficient  (r). 
Upper  Bin  (Mist). 
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Fig-  (>•  5.  Persistence  correlation  coefficient,  r.  lor  3-5/am.  upper 
bin  (Misi) 

St?  =  (l-/h<r^  (6-2) 

where  SE  is  the  standard  error  and  it  is  the  standard  devi¬ 
ation  in  the  predicted  values.  The  resulting  standard 
errors  for  time  0.  listed  in  column  3  of  Table  6.6,  may  be 
compared  with  a  measurement  uncertainty  of  ±2°/cT. 
This  uncertainty  corresponds  to  an  uncertainty  in 
In  (3-5/a in)  (for  the  median  value  of  (3-5/a/w)  in 
the  upper  bin)  of  about  22-38%,  which  compares  reason¬ 
ably  with  the  approximate  standard  error  of  33-55%. 

Column  4  of  Table  6  6  lists  the  inverse  of  the  relax¬ 
ation  constant,  a~\  (for  3-5 uni),  converted  to  hours. 
This  time  is  the  interval  over  which  the  correlation 
coefficient  would  be  expected  to  reduce  to  1/e,  or  .37. 
The  time  interval  a'1  is  short  in  the  summer,  lasting 
about  2  hours  Note  that  it  is  somewhat  longer  in  the 
winter,  indicating  that  the  magnitude  of  infrared  extinction 
at  3-5 n>n  was  more  persistent  in  the  winter  than  in  the 
summer 

The  inverse  of  the  relaxation  constant  a  '  may  also 
be  calculated  for  the  visible,  since  both  r (4)  and  r (60) 
values  were  available  The  resulting  a  1  values,  given  in 
column  5  of  Table  6.6.  are  significantly  longer  than  the 
a  1  values  for  the  3-5 fxm  band.  Thus  the  magnitude  of 
the  visible  aerosol  extinction  apparently  persists  longer 
than  the  magnitude  of  3— 5/u m  aerosol  extinction  for  this 
fog  bin.  That  is,  the  visible  extinction  is  more  stable  than 
the  infrared  extinction.  One  conceivable  cause  for  this 
might  be  that  the  large  droplets  are  more  transient  than 
the  small  droplets 

It  is  interesting  to  compute  the  3-5 fxm  persistence 
correlation  coefficient  for  the  time  duration  of  the  average 
fog  episode.  The  approximate  duration  of  the  average 
mist  or  fog  period  was  calculated  by  determining  the  aver¬ 
age  number  of  consecutive  hours  during  which  the  data  fit 
into  the  upper  bin  This  average  time  interval  is  denoted 
i4  m  column  6  of  Table  6  6.  The  expected  correlation 


Table  6.6.  Related  Persistence  Statistics  for  Upper  Bin  (Mist). 
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coefficient  for  this  average  time  interval  (based  on 
Eq.  6.2)  is  listed  in  column  7  of  Table  6.6.  These  correla¬ 
tion  coefficients  are  very  low,  indicating  that  one  cannot 
expect  much  persistence  or  stability  in  the  infrared  aerosol 
extinction  during  the  average  mist  or  fog  period. 

Thus  these  persistence  calculations  show  that  the 
infrared  extinction  is  relatively  unstable  in  the  mist  or  fog. 
It  varies  on  a  faster  time  scale  than  does  the  visible.  Dur¬ 
ing  the  period  of  the  average  mist  episode,  there  is 
sufficient  variation  in  the  IR  extinction  to  result  in  a  low 
self-correlation  between  the  IR  extinction  values.  These 
variations  could  quite  possibly  be  due  to  short  term  varia¬ 
tions  in  particle  size  distribution. 

These  persistence  calculations  bear  particularly  on 
the  question  of  what  is  causing  the  large  range  in  infrared 
to  visible  ratio.  The  important  point  is  that  the  measure¬ 
ments  do  not  show  a  steady  Ik  to  visible  extinction  ratio 
during  a  given  episode,  with  different  extinction  ratios 
during  other  episodes  on  another  day.  (If  that  were  the 
case  one  might  suspect  that  the  episode  with  the  low  ratios 
was  not  actually  mist  or  fog.)  What  the  measurements 
show  is  a  series  of  episodes  during  which  the  visibility  is 


consistently  reduced,  and  the  relative  humidity  con¬ 
sistently  high;  and  within  each  episode,  the  infrared 
extinction  is  varying  over  a  large  range,  on  a  relatively 
short  time  scale. 

It  is  perhaps  helpful  to  give  examples  of  mist 
episodes.  Figure  6-6  shows  two  long  but  otherwise  fairly 
typical  mist  episodes.  Even  though  the  relative  humidity 
was  over  94%  and  the  visible  extinction  over  1  knr',  both 
the  IR  aerosol  extinction  and  the  IR  to  visible  extinction 
ratio  are  quite  variable.  Even  during  periods  when  the 
visible  extinction  was  over  3  km1,  the  IR  to  visible  ratio 
ranged  from  .2  to  1.2  during  one  episode  and  .2  to  1.0 
during  the  other.  These  two  are  isolated  examples  of  the 
IR  variability  within  the  mist  episodes:  the  persistence 
statistics  show  quantitatively  the  extent  to  which  this 
occurred  generally. 

In  summary,  the  variations  in  infrared  aerosol 
extinction  appear  to  occur  on  a  relatively  short  time  scale. 
In  the  mist  bin,  the  infrared  aerosol  extinction  is  high  a 
significant  portion  of  the  time  but  less  often  than 
expected,  and  the  extinctions  do  not  persist  over  long 
periods. 


Fig.  6-6.  Extinction  variations  within  two  mist  episodes  Dev  15  at  1700  to  Dev  16  at  1000.  1977.  and  Feb  1  from  0300  to  1600.  1978 
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A  possible  future  approach  to  the  problem  of 
predicting  the  high  infrared  aerosol  extinctions  is  discussed 
in  Section  7  For  the  present,  it  appears  that  the  ability  to 
predict  the  extinction  coefficient  in  the  upper  bin  is  prob¬ 
ably  limited  by  short  term  variations  in  infrared  extinction, 
and  that  as  a  result  one  might  profitably  use  a  stochastic 
model  for  this  bin.  Although  the  character  of  the  mist  is 
such  iiiai  one  caiinoi  give  an  exaei  vaiue  of  infrared 
extinction  for  each  value  of  a  predictor,  a  stochastic  model 
still  yields  useful  accuracy,  as  will  be  shown  in  the  next 
section. 

6.7  Estimation  of  Infrared 
Extinction  for  Mist  and  Fog  Cases 

This  section  describes  how  to  use  the  statistics  in 
the  upper  (mist)  bin  to  provide  an  estimation  of  infrared 
aerosol  extinction  for  this  location.  As  described  in  the 
previous  section,  it  was  found  that  for  the  upper  bin 
(unlike  the  lower  bin),  the  ratio  aaer/S  was  better  behaved 
than  the  extinction  aae ,  The  distribution  of  the  aaer/S 
ratio  is  not  a  normal  distribution,  so  as  with  the  lower  bin, 
it  is  necessary  to  use  a  plot  or  table  of  cumulative  fre¬ 
quency  of  occurrence.  The  data  from  the  two  summers 
have  been  combined.  Table  6.7  shows  some  cumulative 
frequency  of  occurrence  values  for  each  season  and  filter. 

These  tables  would  be  used  just  as  the  tables  for  the 
lower  bin  One  would  first  use  the  visible  extinction  and 
relative  humidity  to  determine  that  the  upper  bin  is  the 
appropriate  bin  to  use  Then  one  would  multiply  the 
appropriate  value  of  <xa„/S  from  Table  6.7  by  the  meas¬ 
ured  value  of  visual  extinction  to  determine  the  estimated 
infrared  extinction.  (The  appropriate  value  of  aac,/S  to 
use  would  depend  on  the  application;  on  whether  one 
wanted  the  median  value,  for  example,  or  perhaps  the 
value  that  will  not  be  exceeded  95%  of  the  time  ) 

The  expected  error  in  estimation  for  this  location  is 
shown  in  Fig.  6-7.  These  values  are  based  on  the  50% 
value  of  cioet.'S.  In  Fig.  6-7(a&c),  the  range  for  the  x-axis 
is  the  same  as  for  Fig  6-1.  i.e.  from  0  to  the  median  Sum¬ 
mer  '77  value  of  ( a.MOL  +  a«,o)J-  The  median  error2 
values  are  listed  on  the  figures.  ‘  The  errors  are  less  than 
the  molecular  extinction  about  80%  to  95%  of  the  lime, 
depending  on  the  season. 


Table  6.7.  Cumulative  Frequency.  Ratio  of  Infrared  to 
Visible  Aerosol  Extinction  ( aarr/S ).  Upper  Bin  (Mist ) 
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In  the  remaining  5-20%  of  the  cases,  the  error  can 
become  quite  significant,  as  illustrated  in  Fig.  6-7 (b&d). 
These  figures  illustrate  the  same  data  as  Fig.  6-7(a&c) 
plotted  on  a  different  scale. 

As  with  the  lower  bin,  the  estimated  aerosol  extinc¬ 
tion  applies  reasonably  well  to  any  wavelength  within  the 
measurer!  wavebands,  IT-c  estimate  of  total  extinction 
would  be  obtained  by  using  a  mode!  such  as  LOWTRAN 
to  compute  the  molecular  and  water  vapor  extinctions  at 
the  desired  wavelength,  and  adding  these  to  estimated 
aerosol  extinction. 

6.8  Comparison  of  Measured  and 
Model  Mist  and  Fog  Aerosol  Extinctions 

The  measured  data  in  the  upper  bin  may  be  com¬ 
pared  with  the  LOWTRAN5  model  using  the  ratios  of 
infrared  to  visible  aerosol  extinction  aaer/S.  The  meas¬ 
ured  ratios,  listed  in  Table  6.8,  may  be  compared  with  the 
ratios  inherent  in  the  aerosol  portion  of  the  model.  The 
LOWTRAN5  ratios  for  the  99%  relative  humidity  haze 
models  and  the  fog  models  are  listed  in  Table  6.9.  Figure 
6-8  illustrates  the  model  to  measurement  comparison. 

Comparing  Tables  6.8  and  6.9,  one  can  see  that  at 
least  50%  of  the  measured  ratios  for  each  season  fall  lower 
than  any  of  the  model  ratios  except  the  tropospheric.  In 
fact  the  50%  ratios  for  the  high  bin  are  lower  than  the 
rural,  urban,  or  maritime  ratios  for  any  relative  humidity. 
Thus  many  of  ihe  ratios  are  somewhat  lower  than 
predicted.  The  80%  cumulative  frequency  measured  ratios 
from  Table  6.8  compare  well  to  the  99%  RH  values  from 
Table  6.9.  And  the  95%  or  98%  cumulative  frequency 
measured  ratios  compare  well  to  the  fog  model  ratios. 


Table  6.8.  Measured  Ratios  of  Infrared  to  Visible 
Aerosol  Extinction  (a0„/5).  Upper  Bin  (Mist) 
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Table  6.9.  LOWTRAN5  Ratios  of  Infrared 
to  Visible  Aerosol  Extinction  ( aarr/S ). 
Fog  and  99%  Relative  Humidity  Haze  Models 
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92 
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68 

Rad  log  2 
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24 

Adv  Fog  1 
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1  23 
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1  09 

1  24 
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Thus  in  the  uppe-  bin,  which  is  characterized  by- 
visual  extinctions  greater  than  1  km' 1  and  relative  humidi¬ 
ties  over  94%,  some  of  the  measured  ratios  are  typically 
expected  fog  ratios,  but  many  are  t/pical  haze  ratios  or 
lower. 


One  significant  way  in  which  these  data  differ  from 
the  LOWTRAN  model  is  the  manner  in  which  these  ratios 
in  the  upper  bin  vary.  Ba<=ed  on  LOWTRAN,  one  might 
expect  to  chose  a  model  based  on  the  recent  airmass  his¬ 
tory  or  the  current  conditions.  One  would  for  example 
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Fig.  6-8.  Comparison  of  measured  and  model  ratios  of  infrared 
aerosol  extinction  to  sisible  extinction,  upper  bin  (Mist  and  Fog) 


choose  the  urban  99%  model  if  fog  had  not  formed,  or 
one  of  the  fog  models  if  fog  had  formed.  If  this  informa¬ 
tion  is  lacking,  the  LOWTRAN5  recommendation  is  to 
use  99%  aerosol  models  for  5 of  2  to  4  km  radiation  fog 
for  S  of  4  to  10  km  '.  and  advection  fog  for  S>  10  km"'. 
From  this  model,  one  might  expect  the  ratios  to  be  rea¬ 
sonably  constant  for  a  fair  time  period  In  fog,  for  exam¬ 
ple.  one  might  expect  to  use  ratios  above  1  as  long  as  the 
fog  lasted.  In  contrast  to  this  expectation,  the  OPAQUE 
data  show  that  within  a  given  episode  of  reduced  visibility 
and  high  relative  humidity,  the  infrared  to  visible  ratio  is 
quite  variable,  on  a  short  time  scale  Additionally,  the 
high  ratios  are  not  distributed  with  visual  extinction  in  the 
same  way  one  might  expect  from  the  model.  There  are 
many  high  a0(.r/S  ratios  in  the  5=1-2  km'1  range,  as 
well  as  low  ratios  at  higher  values  of  5.  The  variations  in 
Table  6  8  do  not  appear  to  correspond  to  changes  in  condi¬ 
tions  from  non-fog  to  fog,  but  rather  correspond  to  short 
term  variations  within  the  fog. 


These  results  can  explain  one  reason  that  experi¬ 
mentalists  may  obtain  results  significantly  different  from 
those  predicted  by  LOWTRAN.  Basically,  LOWTRAN 
predicts  one  extinction  for  a  given  condition.  These 
OPAQUE  data  show  that  the  model  extinction  is  a  pretty 
good  estimate  of  the  median  condition,  but  that  the  aero¬ 
sol  extinction  varies  significantly  about  that  median  as  a 
result  of  short  term  variations  in  IR  extinction,  which  may 
be  associated  with  short  term  variations  in  particie  size  dis¬ 
tribution. 

The  upper  bin  data  may  also  be  compared  with  the 
aerosol  portion  of  the  Huschke  model.  For  the  upper  bin, 
the  measured  relative  humidities  are  between  94%  and 
.'00%.  This  corresponds  to  a  range  of  .1  to  1.0  in 
Huschke’s  8—12 /xm  infrared  to  visible  ratio.  As  may  be 
seen  in  Table  6.8,  some  of  the  observed  8-12 fim  ratios 
fall  within  this  range,  but  most  are  lower.  Again,  the 
character  of  the  variations  is  somewhat  different  from  the 
Huschke  model.  Even  in  cases  where  the  relative  humi¬ 
dity  is  constant  for  a  number  of  hours,  and  we  believe  the 
air  to  be  saturated,  i.e.  at  about  100%  relative  humidity, 
the  infrared  to  visible  ratio  is  quite  variable  due  to  the 
short  term  variations  discussed  earlier.  This  model,  as 
well  as  the  LOWTRAN  model,  might  be  improved  by 
using  a  distribution  of  expected  ratios  rather  than  using 
one  expected  ratio. 

6.9  Summary  for  Upper  Bin  (Mist)  Cases 

If  the  visible  extinction  is  greater  than  1  km-'  and 
the  relative  humidity  is  greater  than  or  equal  to  94%,  the 
data  for  this  location  may  be  classed  in  an  upper  bin,  or 
mist  and  fog  bin.  Even  ip  this  bin,  the  infrared  aerosol 
extinction  remains  lower  than  the  molecular  and  water 
vapor  extinction  (about  .5  km-1  for  3-5 fim  and  .25  km-1 
for  8 —  1 2/i./tt)  about  80%  of  the  time.  The  infrared  aerosol 
extinctions  exceed  1  km-1  about  10%  of  the  time.  The 
variations  in  infrared  aerosol  extinction  occur  on  a  rela¬ 
tively  short  time  scale,  so  that  the  self  correlation  in  the 
infrared  extinction  over  the  period  of  the  average  fog 
episode  is  quite  low. 

For  this  location,  the  infrared  aerosol  extinctions  for 
the  mist  bin  may  be  estimated  reasonably  well  about  80% 
of  the  time  through  the  use  of  cumulative  frequency  dis¬ 
tributions  of  the  infrared  to  visible  aerosol  extinction  ratio. 
However,  about  20%  of  the  time  this  method  will  result  in 
large  estimation  errors.  It  should  be  operationally  impor¬ 
tant  to  develop  mechanisms  to  predict  or  detect  these  20% 
of  the  cases  in  which  the  aerosol  extinction  becomes  large. 

7.  SUGGESTED  FURTHER  WORK 
TO  ENHANCE  THE  ANALYSIS 

7.1  Analysis  of  the  OPAQUE  Data 

The  results  presented  in  this  report  are  based  on  a 
large  number  of  data  samples-roughly  7000  hours  worth 
of  data-yet  they  represent  just  one  station  and  three  sea¬ 
sons  (2  summer  seasons  and  a  winter  season).  A  spot 
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check  of  a  few  othar  stations  might  yield  further 
verification  of  the  results  of  this  study.  It  is  reasonable  to 
expect  that  the  basic  technique  of  dividing  the  data  into 
two  bins  should  be  effective  at  other  locations,  however 
one  would  expect  that  the  cumulative  frequencies  of  the 
bins  and  ihe  most  effective  bin  thresholds  are  to  some 
extent  a  function  of  season  and  location. 

It  would  be  particularly  useful  to  look  at  the  distri¬ 
bution  of  aerosol  extinction  in  the  upper  bin  and  to  deter¬ 
mine  the  persistence  of  the  aerosol  extinctions  in  the 
upper  bin  for  other  locations.  Perhaps  the  most  important 
result  of  this  study  was  finding  that  high  infrared  extinc 
tions  are  less  stable  than  expected  in  mist  and  fog  The 
OPAQUE  data  are  ideal  for  determining  the  extent  to 
which  this  is  true  at  other  stations  The  Fall  Netherlands 
data  might  be  particularly  interesting,  since  'an.ssen  and 
\an  Schie  (1981)  note  that  high  extinctions  are  much 
more  frequent  in  the  fall. 

7.2  Furthei  Experimental  Work 

One  of  the  primary  conclusions  of  this  analysis  is 
that  high  infrared  extinctions  occur  less  frequently  than 
expected  in  fog,  and  have  low  persistence  The  variations 
in  infrared  extinction  appear  to  be  the  result  of  short  term 
variations  with  time  and  space.  As  a  result,  if  one  is  to 
predict  the  extinctions  with  greater  accuracy  than  can  be 
obtained  with  a  stochastic  model,  one  scheme  worth  con¬ 
sidering  is  operationally  feasible  techniques  to  detect  the 
presence  of  large  droplets  at  the  exact  time  and  place  the 
prediction  is  desired.  Two  parameters  which  may  be 
strongly  affected  by  the  presence  of  large  droplets  are  the 
wavelength  dependence  of  aerosol  scattering  in  the  visible,, 
and  the  shape  of  the  phase  l unction  for  directional  scatter¬ 
ing  in  the  visible.  Measurements  which  are  affected  by 
these  parameters,  such  as  measurements  of  the  blue/red 
scattering  ratio  or  measurements  of  forward/back  scatter¬ 
ing  ratio  might  provide  means  of  detecting  those  fog  con¬ 
ditions  which  lead  to  high  infrared  aerosol  extinctions.  A 
feasibility  experiment  to  determine  the  efficacy  of  these 
measurements  might  be  well  aavised. 

The  above  experiment  might  provide  means  of 
predicting  the  variations  within  the  upper  bin.  If  it 
becomes  important  tc  predict  the  variations  in  the  lower 
bin,  then  it  would  be  advisable  to  make  additional  meas¬ 
urements  with  much  longer  path  length  instruments  both 
in  the  visible  and  the  infrared  The  lower  bin  suffers  from 
the  instrumental  noise  in  both  the  visible  and  infrared,  but 
longer  path  length  transmissometers  would  measure  this 
region  more  accurately  Since  the  lower  bin  extinctions 
are  low  relative  to  the  molecular  extinction,  however,  it 
should  not  be  necessary  to  resolve  this  variation  in  the 
lower  bin,  for  most  operational  purposes. 

8.  CONCLUSION 

This  report  discusses  an  analysis  of  infrared,  and 
visible  extinctions  recorded  in  the  Netherlands  over  three 
3-month  seasons.  The  data,  recorded  at  hourly  intervals 
during  this  period,  were  gathered  as  a  part  of  the 


OPAQUE  program.  They  are  preliminary  releases  data. 
(See  Appendix  B.) 

It  was  found  that  the  infrared  aerosol  extinctions 
could  effectively  be  sorted  into  a  clear  to  haze  category 
and  a  mist  to  fog  category.  Two  observations  were 
required  to  determine  the  appropriate  category:  the  visual 
extinction  and  the  relative  humidity  The  use  of  either 
parameter  alone  yielded  a  much  less  effective  categoriza¬ 
tion. 

The  haze  category  consisted  of  those  points  with 
visual  extinction  less  than  or  equal  to  1  '.nr1  or  relative 
humidity  less  th^n  94%.  (See  Appendix  B.)  In  this 
category  the  ae_osol  extinction  coefficients  were  con¬ 
sistently  low  with  respect  to  the  water  vapor  and  air 
molecule  extinction  The  median  aerosol  extinctions  were 
about  01  to  02  km  '.  The  extinctions  may  be 
represented  by  a  cumulative  frequency  distribution  of 
aerosol  extinction,  with  a  resulting  error  of  estimate  which 
is  consistently  low.  The  infrared  aerosol  extinctions  in 
this  non-fog  category  are  reasonably  close  to  model  extinc¬ 
tions. 

The  mist  to  fog  category  consisted  of  those  points 
with  visual  extinction  greater  than  1  km  1  and  relative 
humidity  greater  than  or  equal  to  94%.  In  this  category, 
the  aerosol  extinctions  are  typically  about  .1  km  with  a 
significant  number  of  cases  between  1  and  10  km  '  or 
higher.  Thus  the  mist  bin  includes  the  high  values 
expected  in  fog,  however  these  high  values  ate  surpris¬ 
ingly  infrequent  anu  transitory.  The  aerosol  extinction  is 
greater  than  the  water  vapor  and  a>r  molecule  extinction 
(aoout  .5  km  1  at  3-5/xm,  about  .25  km  1  at  8-12 tim)  in 
about  20%  of  the  cases,  and  they  persist  for  short  periods 
relative  to  the  visual  extinction  persistence  and  relative  to 
the  duration  of  the  fog  episodes. 

When  the  aerosol  extinctions  in  this  mist  category 
are  high,  they  match  the  fog  models  reasonably  well,  but 
much  of  the  time  they  are  much  lower.  In  contrast  to  the 
models,  these  data  show  the  existence  of  high  IR  to  visi¬ 
ble  extinction  ratios  at  lower  visual  extinction  thresholds, 
and  show  a  great  variabilif  in  infrared  aerosol  extinction 
within  the  fog 

The  variance  in  infrared  aerosol  extinction  within 
the  fog  was  somewhat  related  to  the  variance  in  visual 
extinction.  As  a  result,  the  infrared  extinctions  for  the 
fog  bin  were  best  estimated  through  the  use  of  the  cumu¬ 
lative  frequency  distribution  of  infrared  to  visible  aerosol 
extinction  ratio  aMJS. 

Although  the  bin  thresholds  and  resulting  cumula¬ 
tive  frequency  distributions  are  applicable  only  to  this  loca¬ 
tion,.  the  sorting  scheme  should  be  applicable  to  other 
locations  as  well.  Perhaps  the  most  important  result  of 
this  work  is  the  observation  that  even  when  both  the 
visual  extinction  and  relative  humidity  are  high,  indicating 
mist  or  fog,  the  infrared  aerosol  extinction  remains  quite 
variable  on  a  short  time  scale,  with  high  infrared  aerosol 
extinction  occuring  frequently,,  yet  less  often  than 
expected  Even  if  much  higher  thresholds  for  visual 
extinction  were  used  to  define  the  upper  bin.  the  high 
infrared  aerosol  extinctions  are  still  a  minority. 
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•  As  a  rcsuii  of  this  work,  one  may  sugpest  two  pri¬ 
mary  future  avenues  of  investigation  First,  interrogate 
the  remaining  OPAQUE  data  to  determine  the  distribution 
and  persistence  of  fie  infrared  aerosol  extinctions  in  the 
mist  at  other  locations  Secondly,  begin  an  experimental 
study  to  determine  the  feasibility  of  detecting  large  water 
droplet  formation  and  thus  the  high  infrared  aerosol 
extinctions  using  measurements  such  as  blue/red  scatter¬ 
ing  ratios.  The  experimental  approach  would  be  to 
explore  measurement  techniques  which  are  operationally 
feasible  and  which  are  theoretically  expected  to  be 
effective  in  detecting  large  water  droplet  formation  and 
thus  high  infrared  aerosol  extinctions. 

v 
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.PPENDIX  B 
Analytic  Update 

As  noted  in  the  text  of  the  report,  this  analysis  is 
based  on  preliminary  release  data  issued  by  the  Nether¬ 
lands  More  recently,  an  improved  data  base  (coded 
"issue  5")  has  become  available  The  relative  humidity 
values  during  fog  have  been  readjusted  from  94%  to  about 
99%  in  the  revised  data  base.  The  problem  with  the  rela¬ 
tive  humidities  discussed  on  p.  6  of  this  report  (last  para¬ 
graph)  has  been  attributed  to  incorrect  maintenance,  and 
the  values  should  be  close  to  100%,  as  suspected  This 
change  in  the  data  base  would  lead  to  a  higher  threshold 
value  than  the  value  of  94%  given  on  p.  10  of  this  report, 
but  would  not  otherwise  affect  the  conclusions. 

The  plots  on  p.  8  show  the  error  bar  associated  with 
1 20%  uncertainty  in  scattering  coefficient  (±20%  of  the 
measured  value).  As  discussed  on  p.  4,  this  was  the 
desired  accuracy  to  be  achieved  by  all  stations  (see  renn 
(1978)  )  More  currently,  the  Netherlands  ass  ises  the 
ALG  accuracy  to  be  about  ±10%  (perhaps  ±  .0%  of  full 
scale). 

The  discussion  of  persistence  calculations,  p.  19. 
notes  that  the  auto-u  elation  /  (r )  is  expected  to  decay 
exponentially,  in  the  absence  of  periodic  variations  It 
should  be  noted  that  periodic  variations  do  occur  in  these 
data,  as  shown  in  Janssen  and  van  Schie  (1982).  however, 
the  observed  period  is  much  too  long  (24  hours)  to  affect 
the  analysis  period  (about  I  hour)  in  this  report.  Sys¬ 
tematic  variations  of  shorter  period  may  in  fact  occur;  this 
is  one  of  the  subjects  undergoing  current  investigation. 

It  may  be  noted  that  the  use  of  the  definition  "visi¬ 
bility  =3,'S”  on  p.  10  is  consistent  with  the  term  "visibil¬ 
ity"  used  by  McIntosh  (1963)  in  the  definition  of  mist 
repeated  on  p  10  McIntosh  states  that  "the  World 
Meteorological  Organization  has  recommended  the  use  of 
the  [contrast  threshold)  value  05.  and  this  figure  should 
be  used  for  the  purposes  of  conversion  of  attenuation 
measurements  to  corresponding  visibility.”  The  relation 
V=3/S  is  derived  using  a  contrast  threshold  .05 

APPENDIX  C 

Further  Comments  on  Infrared  Aerosol 
Extinction  Uncertainties 

Sect  5.2  discusses  four  sources  of  uncertainty 
which  had  to  be  considered  in  analyzing  the  data,  the 
measurement  uncertainty,  the  uncertainty  due  to  non¬ 
simultaneity  of  the  data,  the  uncertainty  associated  with 
the  calibration  procedure,  and  the  uncertainty  resulting 
from  the  molecular  extinction.  This  appendix  contains 
additional  remarks  on  the  interpretation  of  these  uncer¬ 
tainties 


The  first  of  the  four  uncertainties,  measurement 
uncertainty,  is  due  to  detector  noise,  turbulence,  and 
other  such  measurement  problems.  Detector  noise  and 
turbulence  may  be  expected  to  yield  random  errors,  how¬ 
ever  systematic  errors,  both  known  and  unknown,  may 
exist.  The  magnitude  of  the  total  uncertainty  is  difficult  to 
assess.  The  goal  was  to  achieve  ±2%.  Some  error 
sources  tend  to  result  in  a  constant  A  T ,  and  others  tend  to 
result  in  a  constant  A  T/T.  The  actual  error  should  be  a 
fairly  complicated  function.  See,  for  example,  Douglas  et 
al.  (1977).*  We  have  used  A  7"  =  2%  in  order  to  generate 
the  error  bars  in  Fig.  5-4.  Note  that  since  T=e~"r,  we 
have  ( dT/da)=-rT ,  so 

Aa  _  1  A7~ 

a  ra  T 

Thus  if  we  assume  constant  A  T,  A«/«  is  smallest  mid-scale 
on  the  transmissometer.  The  error  bars  for  A«/<*  become 
large  for  T  near  100%  because  «  i»  small,  and  become 
large  for  T  near  0%,  because  T  is  small,  as  shown  in  Fig. 
5-4.  On  the  other  hand,  if  we  assume  constant  A  r/7\ 
then  A  a/a  is  large  only  for  small  a.  Thus,  if  we  had  used 
Af/r=  2%  rather  than  A  7"  =  2%,  the  error  bars  at  low  ,« 
in  Fig.  5-4  would  be  essentially  unchanged,  but  the  error 
bars  at  high  a  would  be  smaller. 

It  should  also  be  noted  that  Table  5.1,  "Aerosol 
Extinction  Uncertainty  for  ±2%  Uncertainty  in  Total 
Transmittance",  does  not  imply  that  we  attribute  the  meas¬ 
urement  error  to  variations  in  aerosol  extinction  What  it 
shows  is,  given  an  error  in  measured  transmittance  due  to 
noise  and  other  sources,-  there  is  a  resulting  error  in  the 


aerosol  extinction  computed  from  that  transmittance.  The 
error  in  the  computed  extinction  is  illustrated  in  Fig.  5-4. 

The  second  of  the  four  uncertainties,  the  uncer¬ 
tainty  due  to  non-simultaneity  of  the  data,  is  not  a  meas¬ 
urement  error.  It  is  simply  a  complicating  factor  in  the 
interpretation  of  plots  such  as  Fig.  5-4.  The  fact  that  the 
measurements  are  not  exactly  simultaneous  adds  some 
uncertainty  to  our  ability  to  interpret  these  plots.  Section 
5.2  discusses  the  magnitude  of  the  time  instability  which 
causes  this. 

The  third  of  the  four  uncertainties,  the  uncertainty 
associated  with  the  calibration  procedure,  is  a  systematic 
error  in  the  transmittances  As  explained  in  Section  5.2,  if 
the  calibration  measurements  were  noise  free,  the  calibra¬ 
tion  constant  would  be  exact,  at  least  to  the  extent  that 
Lowtran  matches  truth.  However  there  is  noise  during  the 
calibration  measurement,  so  the  calibration  constant  has  a 
resulting  small  uncertainty  This  uncertainty  in  the  cali¬ 
bration  constant  results  in  a  systematic  error  in  the 
transmittances.  This  uncertainty  should  be  small  com¬ 
pared  to  the  ±2%  mentioned  before. 

The  fourth  of  the  four  uncertainties  is  the  uncer¬ 
tainty  due  to  the  molecular  extinction  calculation.  In  com¬ 
puting  the  aerosol  extinction,  it  is  necessary  to  extract 
molecular  and  water  vapor  extinction.  If  temperature  and 
relative  humidity  are  not  precisely  known,  the  computed 
molecular  extinction  will  have  an  uncertainty,  and  .here- 
fore  the  computed  aerosol  extinction  will  have  ar  uncer¬ 
tainty.  Note  that  this  is  not  a  measurement  unceriainty  in 
the  infrared  transmittance,  in  fact  the  transmittance  is 
unaffected  by  this  It  is  only  the  final  computed  aerosol 
extinction  which  is  affected. 
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